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PREFACE 

And  God  said,  "Let  there  be  light  . . . and  let  It  be  broadband." 


The  desire  to  Image  distant  objects  Illuminated  by  naturally  occur- 
ring (l.e.,  white)  light  provided  a primary  motivation  for  this  work. 

The  analysis  of  a idilte-llght  shearing  Interferometer  necessitated  the 
development  of  a coherence  model  which  would  adequately  describe  the 
Interference  of  this  so-called  "Incoherent  light."  Analysis  of  the 
Interferometer  Is  presented  as  an  application  of  the  broadband  field 
model  developed  In  Chapter  II. 

I would  like  to  thank  my  sponsor  Donald  W.  Hanson,  Rome  Air 
Development  Center  (OCSE) , for  initially  providing  me  with  this 
Interesting  problem,  and  for  continuing  to  support  me  throughout  the 
study.  I also  gratefully  acknowledge  Professors  Donn  Shankland  and 
Peter  Maybeck  for  serving  as  readers,  and  for  contributing  their  comments 
and  constructive  criticism  to  the  final  draft. 

I would  like  to  especially  thank  Captain  Stanley  R.  Robinson  for 
ably  advising  me  on  the  thesis,  and  for  providing  much  appreciated  gui- 
dance and  motivation  throughout  the  study.  His  Interest  In  the  problem, 
as  well  as  my  personal  development,  had  no  small  part  In  sustaining  my 
enthusiasm  throughout  these  past  months. 


Paul  S.  Idell 
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ABSTRACT 


A free-space  propagation  model  for  broadband  optical  fields  Is 
developed  based  on  a Karhunen-Lo4ve  (KL)  expansion  of  the  time-varying 
portion  of  a coherence  separable  broadband  optical  envelope.  For  long 
characterization  Intervals  It  Is  found  that  the  eigenfunctions  of  the  KL 
expansion  are  approximated  by  complex  exponentials  of  a Fourier  series 
expansion;  the  corresponding  eigenvalues  are  approximated  by  samples  of 
the  temporal  power  spectrum,  sampled  at  the  harmonic  frequencies  of  the 
Fourier  series  expansion.  The  resulting  modal  expansion  provides  an 
Intuitively  simple  Interpretation  of  the  propagation  of  broadband  fields 
and  allows  the  output  field  correlation  to  be  easily  calculated. 

The  propagation  model  Is  applied  to  the  analysis  of  a lateral 
shear  AC  Interferometer,  which  has  been  used  to  measure  the  spatial  phase 
variations  of  a white-light  optical  field  envelope  located  at  Its  Input 
aperture.  It  Is  found  that  the  Interferometer's  operation  for  broadband 
fields  causes  the  phase  which  Is  measured  by  the  Interferometer  to  be 
related  to  a spatially  filtered  version  of  the  aperture  field.  The  effect 
of  this  spatial  filtering  on  the  phase  measurement  Is  studied  for  aperture 
fields  with  arbitrary  spatial  coherence. 

Finally,  the  propagation  model  Is  applied  to  the  shearing  Interfero- 
meter's operation  as  a wavefront  sensor  In  a phase-compensated  imaging 
system,  where  the  phase  of  the  aperture  field  envelope  has  been  disrupted 
by  atmospheric  turbulence.  The  turbulence-induced  phase  perturbation  Is 
modeled  as  a unit-modulus  phase  screen  Introduced  at  the  aperture  plane 
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of  the  Interferometer  for  each  temporal  mode  of  the  field  expansion.  The 
space-time  separability  of  the  source  field  allows  the  visibility  effects 
of  extended  source  distributions  on  the  interferometer's  operation  to  be 
analyzed  separately  from  the  spatial  filtering  effects  of  the  measurement 
process  due  to  broadband  source  emission.  The  interference  fringe  visi- 
bility is  discussed  for  uniform  and  complex  source  radiance  distributions, 
and  the  Interferometer  phase  measurement  is  evaluated  for  typical  phase- 
compensated  imaging  applications. 
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I.  Introduction 

Atmospheric  turbulence,  caused  by  random  fluctuations  in  the  refrac- 
tive index  of  air,  degrades  the  imaging  capability  of  ground-based  astro- 
nomical telescopes.  As  light  from  the  source  propagates  down  through  the 
atmsophere,  turbulence  tends  to  distort  the  shape  of  the  optical  wave- 
front as  well  as  cause  intensity  variations  across  the  wavefront  (Ref 
1:46-48).  Intensity  variations  are  caused  by  a random  lensing  action 
and  give  rise  to  scintillation  and  the  twinkling  of  stars.  The  major 
effect  on  image  quality,  however,  is  due  to  the  random  fluctuations  of 
the  optical  phasefront.  In  long-exposure  telescopic  photography  the 
attainable  resolution  is  turbulence  limited,  rather  than  diffraction 
limited,  to  about  two  seconds  of  arc  (Refs  2 and  3). 

Recently  much  Interest  has  developed  in  the  use  of  predetection 
phase-compensation  for  Improving  the  quality  of  images  distorted  by 
atmospheric  turbulence  (Refs  4,  5,  6,  and  7).  In  the  system  described 
by  Hardy,  ^ (Ref  7) , real-time  phase-compensation  is  accomplished 

by  correcting  the  optical  wavefront  at  the  imaging  system's  input  aperture 
with  a monolithic  piezoelectric  mirror.  The  wavefront  deformation  at  the 
input  aperture  is  determined  with  the  use  of  a lateral  shear  AC  interfero- 
meter (Ref  8)  which  measures  the  wavefront  tilt  of  the  field  at  several 
locations  in  the  aperture  plane.  The  operation  of  this  Interferometer 
as  a wavefront  sensor  in  atmospheric  correction  systems  is  also  discussed 
in  Refs  9 and  10. 
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In  addition  to  atmospheric  compensation  for  Imaging  applications, 
there  exist  other  active  optic  systems  for  which  real-time  correction  of 
wavefront  aberrations  Is  required  (Refs  11,  12,  and  13).  As  a part  of 
these  systems,  the  lateral  shear  AC  Interferometer  may  be  used  to  mea- 
sure wavefront  aberrations  Introduced  by  either  surface  deformations  In 
the  optics  or  fluctuations  In  the  transmission  medium.  A broad  overview 
of  current  active  optic  systems.  Including  descriptions  of  various  wave- 
front  sensors  and  wavefront  correction  devices,  can  be  found  In  Ref  10. 
Several  detailed  articles  on  the  current  theory  and  application  of  active 
optic  systems  also  appear  In  the  March  1977  Issue  of  the  Journal  of  the 
Optical  Society  of  America. 

At  this  time  the  operation  of  the  shearing  Interferometer  men- 
tioned above  has  been  analyzed  with  respect  to  Its  compensated  Imaging 
applications  and  for  limited  classes  of  Input  fields.  A description  of 
the  Interferometer's  operation  for  more  general  aperture  fields  Is  needed. 
Furthermore,  since  sources  of  primary  Interest  In  astronomical  Imaging 
applications  are  Illuminated  with  naturally  occurlng  light  (e.g.,  sunlight), 
the  wavefront  which  must  be  measured  by  the  Interferometer  Is  the  wave- 
front  of  a white-light  (temporally  broadband)  aperture  field.  The  opera- 
tion of  the  shearing  Interferometer  for  extended,  broadband  sources  is 
not  well  understood. 

Problem 

The  object  of  this  study  Is  to  provide  a clear  and  concise  explana- 
tion of  how  the  lateral  shear  AC  Interferometer  measures  the  wavefront 
slope  of  a white-light  aperture  field.  The  effect  of  broadband  light  on 
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the  Interferometer  phase  measurement  is  to  be  determined.  Of  particular 
Interest  to  Its  application  In  compensated  Imaging  systems,  a unified, 
wave-optics  description  of  the  Interferometer  with  regard  to  broadband, 
extended  sources  Is  to  be  given. 

Approach 

This  paper  presents  a description  of  the  lateral  shear  AC  Interfero- 
meter based  on  a model  developed  for  the  free-space  propagation  of  coherence 
separable,  temporally  broadband  optical  fields.  By  definition,  fields  which 
are  coherence  separable  are  those  whose  space-time  correlation  factors  Into 
a product  of  a spatial  correlation  and  a time  correlation.  Coherence 
separability  allows  the  Input  fields  to  be  simply  represented,  and  lends 
tractabllity  to  the  mathematics.  However,  because  the  spatial  part  of  a 
source  field  represented  In  this  manner  is  independent  of  time,  attention 
Is  restricted  to  sources  which  do  not  move  and  whose  Illumination  (or 
luminance)  Is  time- Invariant.  Sources  must  also  be  assumed  to  be  perfectly 
diffuse  (l.e.,  Lambertian),  so  that  their  emission  spectrum  Is  not  a func- 
tion of  viewing  angle. 

By  temporally  broadband  It  Is  meant  that  the  source  emission  spectrum 
is  broad  with  respect  to  Its  medlam  wavelength— for  the  white-light  appli- 
cations which  are  considered  in  this  paper,  the  optical  fields  are  assumed 
to  range  over  the  entire  visible  range,  nominally  O.A  to  0.7  pm.  The 
specific  functional  form  of  the  temporal  power  spectrum  Is  not  specified, 
although  the  spectrum  Is  assumed  to  be  smoothly  varying  over  the  visible 
range.  The  median  wavelength  of  white-light  Is  assimied  to  be  0.55  pm. 
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In  Chapter  II  a free-space  propagation  model  for  broadband  optical 
fields  based  on  a modal  expansion  of  the  time-varying  portion  of  the  field 
Is  developed.  The  temporal  and  spatial  variations  In  the  optical  field 
are  modeled  as  complex  random  processes,  where  the  bandwidth  of  the  sta- 
tionary time  process  corresponds  to  the  bandwidth  of  the  visible  spectrum. 
Spatial  variations  of  the  Input  field  are  assumed  to  be  Independent  of 
time  so  that  the  space-time  field  correlations  are  coherence  separable. 
Using  the  broadband  propagation  model,  the  output  field  correlation  due  to 
broadband  source  fields  Is  calculated  and  output  field  coherence  Is  spe- 
cialized for  sources  with  special  cases  of  spatial  coherence. 

The  broadband  propagation  model  Is  applied  to  the  shearing  Inter- 
ferometer optics  train  In  Chapter  III,  and  the  Interferometer  output 
signal  for  x-shear  Is  determined.  To  simplify  the  analysis,  diffraction 
effects  due  to  finite  lens  size  are  Ignored,  all  optics  are  assumed  to  be 
perfectly  transmissive  and  aberration-free,  and  the  lenses  are  assumed 
to  be  apochromatlc.  The  phase  which  Is  measured  by  the  Interferometer  Is 
shown  to  be  related  to  the  phase  of  the  spatial  part  of  the  complex  aper- 
ture field  envelope.  The  effects  of  white-light  on  the  aperture  field 
phase  measurement,  as  derived  from  the  Interferometer  output  signal,  are 
shown  for  aperture  fields  with  arbitrary  spatial  coherence. 

Finally,  In  Chapter  VI  the  shearing  Interferometer  phase  measurement 
Is  determined  for  an  aperture  field  due  to  an  extended,  white-light  source 
viewed  through  atmospheric  turbulence.  The  source  field  Is  assumed  to  be 
coherence  separable,  corresponding  to  that  reflected  by  a perfectly  diffuse, 
stationary  object  Illuminated  by  white  light.  The  source  Is  also  assumed 
to  be  distributed  and  spatially  Incoherent  as  well  as  being  temporally 
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broadband.  The  optical  phase  distortion  caused  by  atmospheric  turbulence 
Is  modeled  as  a multiplicative  phase  factor  Introduced  at  the  aperture 
plane  of  the  Interferometer  for  each  wavelength.  Validity  of  this  model 
Is  discussed  In  detail.  Fundamental  limitations  on  Interference  fringe 
visibility  and  phase  measurement,  due  to  source  radiance  distribution  and 
emission  spectrum,  are  discussed. 


Functional  Notation 

This  section  Introduces  a standardized  set  of  notation  which  Is  used 
extensively  throughout  the  paper. 

Propagation  Geometry.  All  optical  fields  are  defined  In  planes 
consisting  of  points  r^  - of  a right-handed  rectangular  coor- 

dinate system  (x,y,z).  The  propagation  of  fields  Is  assumed  to  be  In  the 
positive  z-dlrectlon.  Planes  and  points  on  them  are  denoted  by  a sub- 
script Indicating  the  plane's  position  on  the  z-axls.  For  example, 
plane  Pj  consists  of  points  ri  ■ (xi*yi)  for  which  z ••  zj  . Fields 
defined  In  planes  are  denoted  by  the  subscript  Identifying  the  plane, 
e.g.,  optical  field  Ui(ri,t)  Is  located  at  plane  P^  containing  points 
ri  - (xi.yi)  . 

Summation.  Unless  specifically  noted  In  the  text,  all  summations 
range  over  all  Integers  denoted  by  their  Index.  For  example: 


Integration  over  a plane  containing  points  r^  » be  abbre- 

viated by  a single  integral  sign  followed  by  the  differential  dr^  “ dx^dy^ 
Unless  the  llMts  of  Integration  are  explicitly  written,  the  range  of 
Integration  is  doubly-lnfinite.  For  example 


— — f f 

drig(r)  ■ dxi  dyig(xi,yi) 


(1-2) 


denotes  the  doubly-lnfinite,  two-dimensional  Integral  of  the  function  g 
over  points  rj  ■ (xi»yi)  lu  plane  Pj. 

Fourier  Transforms.  Two-dimensional  Fourier  transforms  are  doubly- 
lnfinite  Integrals  defined  in  the  following  way: 


F {g(r)}  - dr  g(r)exp[-j2Tr(f  x+f  y)] 
xy  • j X y 


(1-3) 


where  F {•}  denotes  the  Fourier  transform  with  respect  to  the  variables 
xy  ^ 

X and  y , where  r « (x,y)  ; and  f and  f are  the  transform  varia- 

X y 

bles  associated  with  x and  y , respectively.  The  Inverse  Fourier 
transform  is  defined  in  a similar  manner: 


Ff^f{g(l))  • df  g(f)exp[+j2iT(f  xff  y)] 

XV  X y 


(1-A) 


where  F^  ^ r , denotes  the  Inverse  Fourier  transform  with  respect  to 
X y^  * 

variables  f and  f , when  f ■ (f  ,f  ) : and  x and  y are  the 

X y X y ^ 

space  variables  associated  with  the  transform  variables  f and  f , 

X y 

respectively.  One-dimensional  Fourier  and  Inverse  Fourier  transforms  are 
defined  in  a manner  similar  to  that  above: 


F^{g(x)}  - jdx  g(x)exp[-j2nfx] 
F£tg(f))  ■ jdf  g(f)exp[+j2xxf] 


(1-5) 


(1-6) 
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II.  COHERENCE  PROPERTIES  OF  BROADBAND  OPTICAL  FIELDS 
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In  this  chapter  the  spatial  coherence  properties  of  broadband  optical 
fields  are  Investigated  based  on  a modal  expansion  of  the  time-varying 
part  of  the  field.  Although  results  presented  here  are  derived  for  spe- 
cific application  to  white-light  optical  systems  In  the  visible  (0.4  to 
0.7  pm),  the  models  developed  are  sufficiently  general  and  may  be  applied 
to  a wide  class  of  broadband  optical  and  Infra-red  systems.  The  Input 
fields  are  assumed  to  be  coherence  separable,  and  the  power  spectra  of  the 
time-fluctuations  are  assumed  to  be  smooth  over  all  frequency  ranges  on 
the  order  of  T \ where  T Is  the  characterization  Interval  of  the 
temporal  process. 

Complex  Field  Envelope 

Let  Ui(ri,t)  represent  the  complex  field  envelope  of  a temporally 
broadband  (white-light)  optical  field  propagating  along  the  positive 
z-axls.  The  field  Is  defined  In  a plane  Pi  located  at  z • zi  con- 
taining points  ri  ••  (xi,yi)  . If  the  temporal  fluctuations  of  the  field 
are  centered  at  optical  frequency  fo  , the  scalar  electric  field  fluc- 
tuations of  the  electro-magnetic  field  may  be  written  (Refs  14;  494-499; 
15:12) 


t 
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Ei(ri,t)  ■ Re{Ui(ri,t)exp[-j2TTfot]}  (2-1) 

where  Ex(ri,t)  Is  the  normalized  scalar  electric  field  given  In  volts 
per  meter  per  /ohm  ; Re  {•}  means  "the  real  part  of;”  and  fo  ■ c/Xq  , 
when  c Is  the  speed  of  light  In  vacuum,  and  Xq  Is  the  median  wave- 
length of  the  light. 
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For  simplicity  of  the  model,  the  spatial  variations  of  UjCrj.t) 
are  assumed  to  be  Independent  of  time,  and  the  field  Is  defined  to  be 
coherence  separable; 

Ui(ri,t)  - ui(ri)w(t)  (2-2) 

where  ui(ri)  represents  the  spatial  part  of  the  Input  field  and  w(t) 
Is  the  broadband,  time-varying  portion  of  the  field.  Strictly  speaking, 
this  coherence  model  Is  applicable  to  source  fields  whose  spatial  varia- 
tions are  constant  for  all  time.  This  restriction  can  be  loosened  some- 
what for  the  case  where  the  spatial  variations  are  constant  over  time 
Intervals  comparable  to  the  characterization  time  T of  the  temporal 
process.  As  a convention  ui(ri)  Is  assigned  units  volts  per  meter 
per  i^ohm  ; w(t)  Is  dimensionless. 

The  spatial  and  temporal  parts  of  the  optical  field  ui(ri)  and 
w(t)  are  assumed  to  be  complex,  zero-mean  random  processes  which  may  be 
represented  In  terms  of  their  real  and  Imaginary  quadratures  as 

ui(ri)  - uij^(ri)  + juij.(ri)  (2-3) 

and 

w(t)  - + jwj(t)  (2-4) 

where  uij^(ri)  and  uij(ri)  are  the  real  and  Imaginary  parts,  respec- 
tively, of  the  complex  random  process  ui(ri)  ; and  w (t)  and  w_(t) 

K X 


are  the  real  and  Imaginary  parts  of  w(t)  . If  ui(ri)  Is  assumed  to  have 
Identically  distributed,  but  uncorrelated  quadratures,  the  second  moments 
of  u^(r]:)  may  be  written 


and 


< 


(2-6) 


where  RjCrjjrj)  Is  called  the  spatial  correlation  function,  denotes 

an  ensemble  average,  and  * denotes  the  complex  conjugate.  Similarly, 
if  w (t)  and  w (t)  are  assumed  to  be  identically  distributed  and 
uncorrelated,  the  second  moments  of  the  complex  random  process  w(t)  are 


(2-7) 

(2-8) 


Since  the  exact  form  of  the  spatial  or  temporal  parts  of  the  input 
field  are  not  known,  it  is  reasonable  to  model  them  as  complex  random 
processes — if  for  no  other  reason  than  a lack  of  better  Information.  Also, 
since  for  any  naturally  occurring  phenomenon  it  would  seem  unreasonable 
to  assume  that  one  quadrature  should  take  precldence  over  the  other,  the 
random  processes  are  assumed  to  have  identically  distributed  quadratures. 

Furthermore,  if  w(t)  is  assumed  to  be  stationary,  the  correlation 
of  the  temporal  fluctuations  of  the  field  may  be  written 


R.^(t,t')  - R„(t-t')  - R^(t)  (2-9) 

R^(t)  - F‘^{S^(f)}  (2-10) 

where  t ■ t - t , F^^  {•}  denotes  the  inverse  Fourier  transform  with 
respect  to  f , and  S (f)  is  twice  the  power  spectrum  of  either  w (t) 
or  w^(t)  . Henceforth,  the  temporal  correlation  is  assumed  to  be  sta- 
tionary as  defined  by  Eq  (2-9) . The  space-time  correlation  or  mutual 
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correlation  function  of  the  broadband  field 


Ul(ri,t) 


Is  defined  as 


ri(ri,rl,t,t')  - ^i(ri,t)ut(rf,t'y 

(2-11) 

“ ^l(ri)u*(ri^  ^w(t)w*(t'^ 

(2-12) 

- Rl(ri,ri)R  (t) 
w 

(2-13) 

Therefore,  the  space-time  correlation  of  a coherence  separable  field 
factors  as  the  product  of  the  spatial  correlation  function  and  the  temporal 
correlation  of  the  field. 

For  typical  white-light  applications,  the  spectral  content  of  w(t) 
may  be  considered  to  range  between  0.4  and  0.7  vim,  corresponding  to  a 
temporal  bandwidth  Af  equal  to  3 x 10^**  Hertz.  For  convenience,  the 
power  spectrum  will  be  normalized,  so  that 

JdfS^(f)  - 1 (2-14) 


Thus,  all  the  power  in  the  complex  field  is  assigned  to  the  spatial  part 
of  the  field  envelope  ui(ri)  . A representative  broadband  power  spectrum 
S (f)  and  its  corresponding  correlation  function  R (t)  are  sketched  In 

w W 

Figure  1,  where  B Is  the  nominal  temporal  bandwidth  defined  by 

B-ifif  (2-15) 

the  coherence  time  of  the  complex  temporal  process  Is  approximately  given 
1 


Figure  1.  (a)  A Typical  Broadband  Temporal  Power  Spectrum  S (f)  and 

(b)  Its  Corresponding  Correlation  Function  R (t) 

V 


Modal  Expansion  of  Complex  Field  Envelope 

In  this  section  a modal  expansion  of  the  time-varying  portion  of  the 
complex  field  envelope  w(t)  is  developed.  A Karhunen-Lo6ve  (KL)  expan- 
sion is  used  so  that  the  coefficients  of  the  expansion  are  uncorrelated 
over  the  characterization  interval.  For  long  characterization  times,  it 
is  shown  that  the  basis  functions  of  the  KL  expansion  can  be  approximated 
by  complex  exponentials,  providing  an  intuitively  simple  interpretation 
of  the  propagation  of  broadband  fields. 


m 


Modal  Expansion  of  w(t)  . Consider  the  following  modal  expansion 
of  the  complex  random  process  w(t)  along  a complete  orthonormal (CON) 

fit  fp 

set  of  basis  functions  {<|)^(t)}  over  a finite  time  Interval  ^ 


w(t) 


1.1. m. 

N-+0O 


I w <(>  (t) 

^ n n 
n--N 


T T 
for  - "2  < t <-j 


(2-16) 


T 

2 


where 


w^  ■ Y I w(t)4>*(t) 


(2-17) 


-T 

2 


"l.l.m."  denotes  limit  In  the  mean,  implying  a mean-square  convergence 
of  the  sum  (2-16)  , and  n is  the  Integer  Index  of  the  n*"^  temporal  mode. 
Note  that  {(f^(t)}  Is  a set  of  complex  functions  yet  to  be  specified. 
Also,  since  w(t)  is  assumed  to  be  a zero-mean,  complex  random  process 
{w^}  are  zero-mean , complex  random  variables. 

By  proper  selection  of  basis  functions  {(Ji^(t)  } ' it  is  possible  to 
expand  w(t)  so  that  the  coefficients  of  the  expansion  pair- 


wise uncorrelated: 


^n<')  ” {"n)  (’'n')  “ ° 


(2-18) 


for  n n"  . The  second  equality  comes  about  because  the  coefficients 
are  zero-mean  random  variables. 

A necessary  and  sufficient  condition  for  the  {w^}  to  be  uncorre- 
lated is  that  the  basis  functions  {()>  (t) } are  the  solutions  to  the 

n 

Fredholm  equation  (Ref  16:180): 
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A 


1 

2 


for 


T T 
- 2 < < 2 


J dt'  R„(t.t') 

-T 

2 


(2-19) 


where  {y^}  are  the  real  eigenvalues  associated  with  the  eigenfunctions 
{<Jn(t)}  for  all  Integers  n ; -<n  < n < <>□  , The  expansion  of  w(t)  on  a 

CON  set  of  eigenfunctions  over  an  Interval  yielding  uncorrelated  coeffi- 
cient Is  known  as  a Karhunen-Lo€ve  (KL)  expansion.  If  the  basis  functions 
of  (2-16)  are  solutions  to  (2-19)  then  the  modal  expansion  of 
w(t)  Is  such  an  expansion. 

Results  from  linear  Integral  equation  theory  (Refs  17:122-140; 
18:242-246)  state  that  any  square  Integrable  kernel  R^(t,t')  of 
(2-19)  may  be  expanded  In  a series 


R (t.t')  - y Y (t)  (t>*(t') 
w ^ n n ’^n 

n 

for  -f<t,t'<y  (2-20) 

T T 

where  the  convergence  Is  uniform  for  - ■2<t,t  ^ '2  ' (2-20)  Is 

called  Mercer's  theorem.  It  can  be  shown  (Ref  15:409)  that  If  the 
correlation  function  of  a zero-mean,  complex  random  process  w(t)  can  be 
expanded  In  a form  (2-20) , the  modal  expansion  given  In  (2-16)  will  con- 
verge In  mean-square.  Thus,  for  any  correlation  R (t,t')  which  Is 

w 

T T 

continuous  and  bounded  for  - '^<t,t'*<-j  , the  modal  expansion  of  w(t) 
given  In  (2-16)  will  converge  In  mean-square  to  the  process  w(t)  . 

Eigenfunctions  and  Eigenvalues  for  Long  T.  For  stationary  random 
processes  characterized  over  long  time  Intervals  [-  ]»  It  can  be 

shown  (Ref  16:205-207)  that  the  eigenfunctions  {<^^(t)}  and  associated 
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eigenvalues  {y^^}  which  are  solutions  to  the  Fredholm  equation  (2-19) 
can  be  approximated  by 


for 


T ^ T 
2<t<2 


and 


(ti^(t)  = exp[+j2Tr^t] 


Y - i S (-) 
^n  T w'‘  T •' 


(2-21) 

(2-22) 


where  T is  the  characterization  interval  in  seconds  and  S (^)  is  the 

w T 

power  spectrum  of  the  complex  random  process  w(t)  , defined  in  Eq  (2-10) , 
sampled  at  frequencies  Hertz. 


The  magnitude  of  T needed  for  the  validity  of  the  approximation 
depends  on  how  quickly  changes  near  frequency  f " ^ • for 

smooth  spectra,  long  T means  long  compared  to  the  reciprocal  bandwidth  of 
the  fluctuations  of  the  optical  envelope; 

T » ^ (2-23) 

For  white-light  applications  the  bandwidth  B = 1.5  x 10^**  Hertz  . There- 
fore, for  a characterization  time  T much  greater  than  6.67  x 10*'^ 
seconds,  the  {<Ji^(t)}  of  the  expansion  over  time  interval  of  length  T 
become  the  complex  exponentials  of  a Fourier  series  expansion,  and  the 
eigenvalues  corresponding  to  the  eigenfunctions  become  samples  of 

the  power  spectrum  evaluated  at  the  harmonic  frequencies  of  the 

Fourier  series  expansion. 

Thus,  for  long  characterization  time  T , the  eigenfunctions  of  the 
KL  (modal)  expansion  (2-16)  may  be  approximated  by  complex  exponentials. 


and  the  expansion  for  w(t)  becomes 


T ^ T 
tor  ~ "2  ^ ^ 


where 


w(t)  • w e 
n 


T 

2 


+j2TrYt 


1 

w = ^ 
n T 


dt  w(t)e 


-jZTT^t 


(2-24) 


(2-25) 


-T 

2 

and  the  statistics  of  the  process  w(t)  are  such  that  the  sum  is  assumed 

to  converge  in  mean-square  as  discussed  earlier. 

With  regard  to  the  above  approximation,  it  may  be  noted  that  if  the 

random  process  is  expanded  in  a Fourier  series,  it  can  be  shown  (Ref  19:94) 

that  the  coefficients  of  the  expansion  {w^}  become  uncorrelated  as  the 

expansion  Interval  T gets  long. 

Eigenvalue  Statistics.  The  expected  value  of  the  energy  of  w(t) 

X X T 

in  time  interval  defined 

T 
2 


dt  w(t)w*(t^ 


-T 

2 


I 


(2-26) 


(2-27) 


where  the  modal  expansion  for  w(t)  has  been  used.  Using  Mercer's 
theorem  (2-20) , the  mean  energy  of  the  process  for  long  characterization 
time  T is 


( 


I S„(f  ) 


(2-28) 


Equating  each  term  of  sums  (2-27)  and  (2-28)  and  using  Eq  (2-18)  yields 


) 6 


nn 


(2-29) 
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where  6 ^ is  a Kronecker  delta.  Thus,  for  long  characterization  inter- 

nn 

vals,  the  mean-square  value  of  each  expansion  coefficient  seen 

to  be  just  a sample  of  the  temporal  power  spectrum  evaluated  at  the 
coefficient's  harmonic  frequency. 

Modal  Expansion  of  the  Optical  Envelope.  From  Eqs  (2-2)  and  (2-24), 
the  complex  field  envelope  of  a coherence  separable,  temporarily  broad- 
band optical  field  in  a plane  z = zj  can  be  written 

_ _ 

Ui(ri,t)  » I u(ri)  w^e 
n 

for  (2-30) 


where  the  complex  field  envelope  is  taken  with  respect  to  a carrier  at 
optical  frequency  fg,  as  in  (2-1).  Furthermore,  its  mutual  correlation 
function  can  be  written 


ri(ri,ri,t,t') 


^ (ri,t)ut(rf,t")^ 

I I ^ 1 (7i ) ut (r 


exp[j  ^ (nt  - n't')] 


(2-31) 


(2-32) 


T T 

for  " that  the  modal  expansion  of  w(t)  allows  the 

ensemble  average  over  time  sample  functions  in  the  temporal  correlation 
to  be  replaced  by  a sum  of  expected  values  of  random  variables. 

Since  the  modal  expansion  coefficients  have  been  chosen  so  as  to  be 
pairwise  uncorrelated  (Ref  Eq  (2-29)),  Eq  (3-32)  may  be  further  simplified 

ri(ri,ri,t,t')  - Rl(ri,ri)exp[j2iT  Y (t  - t')] 

n 

for  • (2-33) 
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The  n temporal  mode  of  the  optical  envelope  Is  defined  by  the  n 
term  In  the  series  (2-30) 

_ +j2Tr|t 

Uijj(ri,t)  = ui(ri)w^e 

for  (2-34) 


( 


where  Ui^(ri,t)  Is  the  baseband  representation  of  a monochromatic  optical 
field  at  frequency  f^  = (fo  - y)  » ui(ri)  Is  the  wavelength-independent 
spatial  part  of  the  complex  field;  and  (w^)  Is  the  n^^  coefficient  of  the 
modal  expansion  for  w(t)  defined  In  Eq  (2-25) . 

The  Importance  of  the  modal  expansion  for  the  complex  field  envelope 
given  In  Eq  (2-30)  Is  the  following.  Since  Maxwell's  equations  governing 
propagation  of  electromagnetic  waves  for  homogeneous.  Isotropic  media  are 
linear  In  space  and  time,  each  n^''  term  In  the  expansion — or,  n^^  temporal 
mode — propagates  Individually.  In  terms  of  optical  fields  this  means  that 
each  temporal  mode  may  be  propagated  as  a monochromatic  wave  at  frequency 
f^  * ffl-  Y • superposition,  the  output  broadband  field  Is  simply  the 
sum  of  all  the  Individually  propagated  modes. 

The  free-space  propagation  of  broadband  optical  fields  using  the  modal 
expansion  developed  above  Is  discussed  In  the  next  section. 


w 


Free-Space  Propagation  Model  for  Broadband  Optical  Fields 

Consider  two  parallel  planes  and  P2  separated  by  a free-space 

propagation  path  of  length  Z as  shown  In  Figure  2.  The  propagation  of 
monochromatic  fields  from  Input  plane  at  z - to  output  plane 

P2  at  z ■ Z2  la  governed  by  the  Huygens-Fresnel  Integral  (Ref  20:60) 
which  may  be  written 


r - - - - 

“2(^2)  ■ "Jx^  J <*riui(ri)exp[j|j  |r2-ri|2] 


(2-35) 
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Figure  2.  Free-Space  Propagation  Geometry 


where  ui(ri)  Is  the  complex  field  amplitude  of  the  Input  field,  U2(r2) 

Is  the  field  amplitude  of  the  output  field,  ri  ■ (xi»yi)  Is  located  In 

» ^2  “ (*2»y2)  is  located  In  P2,  ko  * » ^0  i®  wavelength 

Aq 

of  the  light  and  Z ■ Z2  - zi  Is  the  distance  from  plane  Pi  to  P2  In 


In  Eq  (2-35)  the  paraxial  (Fresnel)  approximation  has  been  used.  For 
many  applications,  the  optical  fields  of  Interest  are  confined  to  a region 


about  the  z-axls  whose  maximum  linear  dimension  Is  small  compared  to  the 
propagation  distance  Z . For  these  situations,  the  Fresnel  approximation 
to  the  Huygens-Fresnel  Integral  Is  very  good.  This  form  also  allows  the 
computations  associated  with  wave  propagation  and  diffraction  to  be 
greatly  simplified,  allowing  a "systems"  type  representation  of  these 
effects. 

For  a monochromatic  field  of  wavelength  the  Huygens-Fresnel 

Integral  (2-35)  can  be  rewritten 


and  ui^(ri)  and  U2^(r2)  represent  the  complex  Input  and  output  field 

amplitudes  at  wavelength  X , 

n 

Using  the  modal  expansion  for  a temporally  broadband  optical  field 
given  In  Eq  (2-30),  the  output  field  mode  U2^(r2,t)  due  to  a monochromatic 
temporal  mode  at  optical  frequency  f “ * fo  - ^ at  the  Input  plane 

tl  A L 

n 


can  be  written 


^'^•4 


where 


(2-39) 


“ ^0  e'^n) 

«<") 

n 


C(n) 


1 


n 

foT 


(2-41) 


and  Ui^(ri,t)  is  the  n temporal  mode  of  broadband  field  envelope 
Ul(ri,t)  located  at  plane  Pj  given  previously  in  Eq  (2-34). 

Note  that  both  input  and  output  field  modes  are  baseband  representa- 
tions of  the  n^^  temporal  mode  at  the  input  and  output  planes,  respectively, 
taken  with  respect  to  an  optical  carrier  at  center  frequency  fo  ** 

Aq 

For  broadband  (white-light)  fields,  the  output  field  U2(r2,t)  can 
be  written  as  the  sum  of  the  individually  propagated  temporal  modes: 


( U2(r2,t)  - I U2jj(r2,t) 

n 

for  -f<t<y  (2-42) 

where  U2^(r2,t)  is  given  in  Eq  (2-38)  and  the  sum  converges  in  the  mean. 
For  reference,  Eq  (2-42)  is  called  the  free-space  propagation  model  for 
temporally  broadband  optical  fields. 

Field  Coherence  Properties  of  Broadband  Sources 

The  purpose  of  this  section  is  to  calculate  the  correlation  of  an 
optical  field  at  an  output  plane  due  to  a temporally  broadband  source  with 
arbitrary  spatial  coherence.  First,  the  output  field  correlation  is  cal- 
culated using  the  broadband  field  propagation  model  developed  in  the 

C 
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^ previous  section.  Using  statistical  models  for  the  spatial  coherence  of 

the  source,  results  for  the  spatial  correlation  of  the  output  field  are 
derived  for  spatially  coherent  and  incoherent  sources.  Results  for  the 
spatial  correlation  of  broadband  light , derived  using  the  broadband  field 
propagation  model,  are  compared  to  well-known  results  for  quasl-monochromatic 
fields. 

Output  Field  Correlation.  The  source  field,  defined  In  plane  Pj  , 
is  assumed  to  be  coherence  separable  so  that  the  Input  field  envelope 
may  be  written  as  in  Eq  (2-30) : 

_ _ +j2lT^t 

Ui(ri,t)  » I ui(ri)w  e (2-43) 

n “ 

where  ui(ri)  represents  the  spatial  part  of  the  complex  field  and  the 
modal  expansion  for  temporally  broadband  (white-light)  fields  has  been 
used. 

Using  the  broadband  field  propagation  results  given  in  Eq  (2-42),  the 
output  field  at  plane  P2  due  to  a temporally  broadband  optical  field  at 
Pi  is 

U2(r2ft)  ~l  ® 

n n 

— — ^ _ 

/dri  ui(ri)exp[j |r2-ri  I 2]  (2-44) 

where  Z is  the  distance  between  planes  Pi  and  P2  In  meters. 


The  correlation  of  the  broadband  output  field  can  be  written 


r2(r2,r2,t,t')  » 


^U2(r2,t)U*(;£,t'^ 

^ f (t-t')] 


|dr.  ^ 


drf  Ri(ri,ri) 


exp[j^(|r2-ri|2-  |rf-rr|2)] 


(2-45) 


(2-46) 


T T 

for 


I(X^Z)‘2dri  dri  ri^(ri,rf,t-t') 

exp[j  rf  (|72-rip-  Irz-FiP)] 


(2-47) 


where 


rin(ri,ri,T)  - is  ( ^ )»!  (^i  ,ri)exp[j2Tr  ^ t ] (2-48) 


is  the  n temporal  mode  of  the  correlation  of  the  source  field  given  in 
Eq  (2-33).  The  output  field  correlation  r2(r2,r2,t,t ')  or  the  mutual 
correlation  function  expresses  the  mutual  coherence  of  light  fluctuations 
at  points  r2  and  rl  in  the  output  plane,  where  the  fluctuations  ait 
r2  are  taken  at  time  t , and  those  at  rj  are  taken  at  time  t". 

Spatial  Correlation  of  Output  Field.  The  spatial  correlation  of  the  out- 
put field  is  calculated  from  Eq  (2-47)  when  t-t'-T“0: 


^2(^2. 1^2)  ” r2(r2.r2,T“0) 


(2-49) 


- I (A^Z)“^|dri  I dri  ri^(ri,ri) 


( • 

1 .( 


T T 
for  ~ 2^^  *2 


exp[j-57  (lr2-riP- |r2-ri|2)] 


(2-50) 


V. 


where 


TljjCri.ri)  - ri^(ri,ri,T-0)  (2-51) 

- Y Sw(f ) Ri(7i,7f)  (2-52) 

The  output  spatial  correlation  function  r2(r2,r2)  Is  also  called  mutual 
Intensity  (Ref  14:508-510)  of  the  output  field.  It  expresses  the  spatial 
coherence  of  light  at  two  points  r2  and  r^  In  the  output  plane  when 
both  points  are  considered  at  the  same  time.  For  white-light  sources 
Eq  (2-50)  Is  directly  related  to  the  fringe  visibility  of  a Mlchelson 
stellar  Interferometer  (Ref  22:30)  for  extended,  white-light  sources  with 
arbitrary  coherence  properties. 

Let  the  correlation  of  the  spatial  part  of  the  source  field  be  written 

Rl(ri,ri)  - ili(ri,ri)ui(ri)ut(ri)  (2-53) 

where  p](ri,rf)  Is  defined  to  be  the  complex  degree  of  coherence  of  the 
source,  such  that  spatially  coherent  sources  are  defined  as  those  for 

Pl(ri,ri)  - 1 (2-54) 

spatially  Incoherent  sources  are  defined  so  that 

fil(ri.7i)  - 6(ri-ri)  (2-55) 

Using  Eq  (2-53)  for  spatially  coherent  but  temporally  broadband  sources, 

the  spatial  correlation  of  the  output  field  may  be  written  from  Eq  (2-50): 
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^2(r2»r2)|  coherent 


T T 

T < t < ■:r 
2 2 


~ I (.X  S (^) 

^ n ' T w'  T ' 
n 


f — — ^ 

J dri  ui(ri)exp[j-^  lr2-ri|2] 

J drf  ui(ri)exp[-j^  |r2-rf|2] 


(2-56) 


For  broadband, spatially  Incoherent  sources  Eq  (2-50)  becomes 


^2(’^2»r2)  incoherent 
source 


T T 
for  - Y < t 


^ s„(f) 

n i w X 

- 


'n  f — — 

dri  Ii(ri)exp[-j-^  ri*(r2-r2)] 


(2-57) 


where 


'('n  “ 2Z  0^2p-|r2p) 


(2-58) 


Il(ri)  - ui(ri)ut(ri) 


(2-59) 


Is  defined  to  be  tHe  Intensity  of  the  source  in  watts  per  meter  squared. 
Note  that  Eq  (2-57)  Is  an  extension  of  the  Van  Clttert-Zernlcke  theorem 
(Ref  14:508-511)  for  temporally  broadband  sources. 

In  particular.  If  the  source  Is  radiating  quasl-monochromatlc  light 
at  wavelength  Xq  , the  temporal  power  spectrum  may  be  approximated  by 


S„(y)  - T 6(n-0) 


(2-60) 
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and  Eq  (2-57)  yields  the  mutual  Intensity  for  an  extended,  incoherent, 
quasl-monochromatic  source  derived  by  Born  and  Wolf  (Ref  14:509); 


r2(r2.r2) 


quasi-monochromatlc , 
Incoherent  source 


where 


» (XqZ)”^  e^"^®  jdn  Ii(ri)exp[-j  ^ ri»(r2-r2)] 

(2-61) 

(|r2l^-|r2p) 

(2-62) 


and  Xq  Is  the  median  wavelength  of  the  source  radiation. 

Note  that  for  each  temporal  mode  n the  mutual  Intensity  of  the 
light  at  the  output  plane  is  just  as  prescribed  by  the  Van  Cittert-Zernlcke 
theorem.  The  spatial  correlation  for  white-light,  extended  sources  (2-57) 
is,  therefore,  the  sum  of  all  the  mutual  intensities  for  each  temporal 
mode.  Also,  since  the  spatial  variations  of  the  source  field  are  assumed 
to  be  wavelength- Independent,  the  contribution  from  each  temporal  mode  is 
weighted  by  the  temporal  power  spectrum  evaluated  at  the  appropriate 

harmonic  frequency. 

Coherence  Length  for  Broadband  Optical  Fields 

Bringing  the  sum  over  temporal  modes  inside  the  integrals  over  r] 
and  writing 

- to 
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T 


allows  Eq  (2-50)  to  be  rewritten  in  the  following  form: 


r2(r2,r2)  “ (XqZ) 


-2 


dri  jdri 


Rl(ri,ri) 


ex 


p[j  ||  (|r2-  riP-  r2-rf|2)] 


- T ^ T 

for  -2<^<2 


where 


(2-64) 


R,(t)  ” ^ I 5^(n)S^(  Y)exp[-j2irT(^)] 


(2-65) 


For  wide  process  bandwidth  B and  long  characterization  time  T , the  sum 
over  temporal  modes  may  be  approximated  by  an  Integral  over  v (Ref  16:207) 


R (v)  = 

w 


dv  C^(v)S  (v)exp[-j2TrTv] 
w 


(2-66) 


F {C2(v)s  (v)} 


w 


(2-67) 


where  C(v)  ■ (1  - t ) and  v has  been  substituted  for  ^ as  the 

IQ  T 

integration  variable. 

Eq  (2-64)  Indicates  that  the  spatial  correlation  of  a broadband  field 

can  be  calculated  from  a two-fold  spatial  Integral  over  (1)  the  spatial 

correlation  of  the  source  Ri(ri,rr)  and  (2)  the  contribution  from 

differential  path  delays  of  all  the  temporal  modes,  written  as  R (•)  • 

w 

As  shown  In  Eq  (2-67)  the  differential  path  contribution  to  the  spatial 
correlation  can  be  approximated  by  the  Fourier  transform  of  a weighted 
power  spectrum.  The  relationship  of  the  quadratic  term  C^(v)  to  a 
typical  power  spectrum  Is  shown  In  Figure  3. 
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Figure  3.  Relationship  Between  Quadratic  Term  C^(v)  and  the  Temporal 

Power  Spectrum  S (v) 
w 


Expanding 


C2(v) 


(2-68) 


and  using  the  differentiation  theorem  for  Fourier  transforms  (Ref  23:36), 

R^(t)  may  be  written  In  terms  of  the  first  two  derivatives  of  the  temporal 

correlation  function  R (t)  : 

w 


R (t)  - R (t)  - 
w w j2nfo 


(jwjrfoT^  3t2 


R (t)  (2-69) 


R (t)  + - 

W TTfo 


3t  “ (2Trfo)^ 
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3t^ 


R^(t) 


(2-70) 


%rtiere  the  correlation  of  the  stationary,  wideband  process  w(t)  Is  given 
(2-10).  If  the  bandwidth  B of  S (f)  Is  narrow  enough  so  that 


for  all  temporal  modes  with  significant  energy,  then  C^(v)  = 1 , and 

R (t)  = R (t)  (2-72) 

w w 

because  the  correlation  functions  considered  here  are  even 

F {S  (v)}  - F~^{S  (v)}  - R (t)  (2-73) 

V W V w w 

For  any  wideband  power  spectrum  of  practical  Interest  tends 

to  drop  off  in  magnitude  for  increasing  argument  (see  Figure  1).  This 
tendency  will  cause  a "windowing"  effect  in  the  Integral  (2-64) . In 
terms  of  spatial  coherence,  there  exists  regions  in  Pj  (consisting  of 
points  ri  and  v\  ) and  P2  (consisting  of  points  r2  and  rz  ) for 
which  the  spatial  correlation  function  r2(r2,r2)  is  negligible.  If  the 
spatial  correlation  of  the  source  is  stationary  for  all  separations 
|ri-ri|  in  the  source  plane,  the  output  coherence  is  also 

stationary  and  may  be  written  r2(|r2-r2|)  . The  separation  = Ir2“r2| 

for  which  the  output  correlation  function  has  significant  value  is  called 
the  coherence  distance. 

Let  the  spatial  correlation  of  the  source  be  modeled  as  in  Eq  (2-53) 
and  assume  the  complex  degree  of  coherence  of  the  source  Pi  is  stationary 
in  spatial  coordinates  (homogeneous) 

RlC^i.ri)  - Pi(pi)ui(ri)ui(ri)  (2-74) 

where  Pi  ■ Ri-Ff]  (2-75) 
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The  spatial  correlation  of  the  field  U2(r2,t)  can  then  be  written 
r2(’^2»J^2)  = (^qZ)  jdrf  Pi(pi)ui(ri)ui(rf) 

(|r2-Fi|2-|r2-ri|2)] 
exp[j  ^ (|r2-ri |2-|r2-ri |2)] 

for  (2-76) 

Special  cases  of  the  output  spatial  correlation  function  for  source 
fields  having  various  combinations  of  emission  spectra  and  spatial 
coherence  are  given  in  Table  I.  Monochromatic  fields  are  assumed  to 
have  impuJslve  emission  spectra  as  given  in  Eq  (2-60).  The  spatial 
coherence  of  coherent  and  incoherent  sources  are  modeled  by  Eqs  (2-54) 
and  (2-55),  respectively. 

In  summary,  it  must  be  noted  that  the  output  field  correlations 
listed  in  Table  I are  calculated  for  source  fields  which  are  assumed  to 
be  coherence  separable.  Source  fields  which  can  be  represented  in  this 
manner  are  restricted  to  those  emitted  by  a perfectly  diffuse  and  sta- 
tionary object  whose  spatial  characteristics  are  independent  of  time. 
Furthermore,  since  the  output  field  correlation  (Ref  Eq  (2-47))  cannot, 
in  general,  be  separated  into  a product  of  a spatial  correlation  and  a 
temporal  correlation,  the  output  field  due  to  a coherence  separable 
source  field  is  not  coherence  separable. 

( 
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III.  APPLICATION  OF  BROADBAND  FIELD  MODEL  TO  THE 
ANALYSIS  OF  A WHITE-LIGHT  SHEARING  INTERFEROMETER 


A lateral  shear  AC  interferometer  has  been  used  in  the  measurement 
of  the  phase-front  of  an  optical  wave.  This  chapter  presents  a descrip- 
tion of  this  interferometer's  operation  for  broadband  optical  fields  with 
arbitrary  spatial  coherence  based  on  the  propagation  models  for  broadband 
optical  fields  developed  in  Chapter  II.  The  study  of  the  shearing  inter- 
ferometer is  Intended  to  demonstrate  the  use  of  the  broadband  coherence 
model  in  the  analysis  of  a broadband  optical  system.  The  fundamental 
limitations  Imposed  on  the  wavefront  measurement  by  the  broadband  nature 
of  the  optical  field  are  sought.  As  such,  the  analysis  ignores  .all  dif- 
fraction effects  due  to  finite  aperture  size  and  all  lenses  are  assumed  to 
be  perfectly  transmissive  and  aberration-free. 

The  first  section  presents  a general  discussion  of  the  interferometer 
optics,  summarizing  results  reported  by  Refs  7 and  10.  The  sections  that 
follow  present  the  analysis  of  the  Interferometer  using  the  broadband 
coherence  model.  The  phase  which  is  measured  by  the  Interferometer  is 
shown  to  be  related  to  the  phase  of  the  complex  part  of  the  aperture  field 
envelope.  Results  of  this  chapter  show  the  effect  of  white-light  optical 
fields  on  the  Interferometer's  phase  measurement.  The  median  wavelength 
of  white-light  Xq  Is  assumed  to  be  O.SS  pm. 

Description  of  Interferometer  Optics 

The  lateral  shear  AC  Interferometer  has  been  used  to  measure  the  local 
slope  of  an  optical  wavefront  (Refs  7,  8,  9,  and  10).  As  shown  in  Figure  4, 
the  incoming  wave  is  beamspllt  into  two  similar  channels,  each  consisting 
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Figure  4.  Two-Channel  AC  (Heterodyne)  Shearing  Interferometer 


of  a pair  of  lenses,  a rotating  radial  grating,  and  a detector  array.  The 

purpose  of  each  channel  Is  to  measure  the  wavefront  slope  along  one  of  the 

two  coordinate  axes  of  the  measurement  plane.  Output  signals  from  the 

X-  and  y-  channels  can  then  be  used  to  reconstruct  a phase  map  of  the  two- 

dimensional  wavefront  sensed  by  the  Interferometer. 

The  optics  In  each  channel  are  configured  as  In  Figure  5.  The  two- 

dimensional  optical  field  In  the  Input  aperture  plane  P Is  denoted 

a 

( 
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Figure  5.  One-Channel  Interferometer  Optics 


U^(r,t)  , where  r^  ■ points  in  plane  . The  input  aper- 

ture of  each  channel  coincides  with  the  measurement  plane  shown  in  Figure  4, 

so  that  the  interferometer  measures  the  phase  of  U (r  ,t)  . 

& A 

The  input  field  is  focused  onto  a rotating  radial  grating  G in 
Plane  P by  a diffraction  limited,  chromatically  compensated  (apochromatlc) 

O 

lens  Li  having  focal  length  Fi  . By  diffraction  limited  and  apochromatlc 
it  is  meant  that  the  lens  will  focus  a plane  wave  located  at  P^  to  a dif- 
fraction limited  spot  at  plane  P for  all  wavelengths.  The  field  at  the 

g 

grating  U^(r^,t)  is  chopped  by  the  moving  grating  and  focused  onto  a 
detector  array  at  plane  P^  by  lens  L2  . Lens  L2  has  focal  length  F2 
and  is  assumed  to  be  apochromatlc  and  diffraction  limited. 

( 
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The  field  chopped  by  the  grating  is  diffracted,  forming  multiple 
Images  of  the  input  field  U (r  ,t)  on  the  detector  plane.  Each  dif- 

d d 

fracted  image  is  shifted,  or  sheared,  with  respect  to  the  next  by  an 
amount 


where  s is  the  shear  distance 

F2  Is  the  focal  length  of  lens  L2 
d is  the  line  spacing  of  grating  (grating  period) 

X is  the  wavelength  of  light 

The  direction  of  shear  depends  on  the  coordinate  direction  the  rotating 
grading  cuts  the  field  in  plane  . 

The  radial  grating  is  a circular  glass  disc  with  alternating  clear 
and  opaque  radial  lines  extending  from  near  the  center  of  the  disc  to  its 
outer  edge.  If  the  grating  has  N opaque  lines,  the  line  spacing  at  a 
distance  R from  the  center  of  the  disc  d ■■  2itR/N  . Since  the  grating 
produces  a square-wave  transmittance  at  any  radius  it  is  also  called  a 
radial  Ronchl  grating.  For  the  x-channel  of  the  interferometer,  the  input 
field  is  focused  onto  the  grating  so  that  the  Ronchl  rulings  move  across 
the  spot  in  the  x-dlrectlon.  The  resulting  diffraction  causes  the  Images 
in  the  detector  plane  to  shear  along  the  x-axls.  In  the  y-channel,  the 
radial  lines  cut  the  spot  along  the  y-axls  producing  a y-shear  in  the 
detector  plane.  Figure  5 gives  a representation  of  the  lateral-sheared 
images  of  the  input  field  for  the  x-  and  y-channels  showing  the  zero,  -1, 
and  -fl  diffracted  orders.  Each  channel  is  constructed  so  that  the  focus 
point  on  the  grating  is  variable  and  the  shear  distances  in  the  x and  y 
directions  may  be  varied  Independently. 
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Figure  6.  Representation  of  Lateral  Sheared  Images  In  Detector  Plane 
for  (a)  x-shear  and  (b)  y-shear  Showing  the  +1,  0,  and  -1  Diffracted 
Orders. 


As  the  gratings  rotate,  the  Intensity  pattern  at  each  point  In  the 
detector  plane  Is  modulated,  producing  an  AC  signal  at  the  output  of  each 


detector.  The  fundamental  modulation  frequency  f^  Is  given  by 

f - M 

s d 2it 


(3-2) 


where 


V Is  the  linear  speed  of  the  radial  grating  at  radius  R In 

O 

centimeters  per  second 

d Is  the  grating  period  at  radius  R In  centimeters 
N Is  the  number  of  opaque  lines  In  the  grating 


6 Is  the  angular  speed  of  the  grating  In  radians  per  second 
As  It  Is  shown  later,  the  fundamental  frequency  corresponds  to  the  Inter- 
ference of  the  zero  and  +1,  with  the  zero  and  -1  diffracted  orders  In  the 


► 

f 
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shaded  region  of  the  detector  plane  shown  In  Figure  6.  If  the  square-wave 

Ronchl  grating  Is  used,  no  other  combination  of  diffraction  orders  produces 

an  Interference  pattern  modulated  at  f^  . It  Is  to  be  shown  how  the 

detector  plane  Interference  pattern  modulated  at  f^  relates  to  the  phase 

of  a broadband  Input  field  U (r  ,t}  . 

a a 


Calculation  of  Detector  Field 

The  propagation  geometry  for  the  x-channel  of  the  lateral  shear  AC 
Interferometer  Is  shoim  In  Figure  7.  The  broadband  field  U (r  ,t) 

3 d 

located  In  the  Input  aperture  (measurement  plane)  P of  the  Interfero- 

a 

meter  can  be  expanded  Into  temporal  modes: 


U (r  ,t)  = y u (r  )w  e 
a a’  ^ an  a n 
n 


+j2iT^t 


T T 
for  - \ 


(3-3) 


where  spatial  part  of  the  aperture  field  and  the  sum  over 

temporal  modes  n converges  In  the  mean-square.  Note  that  the  aperture 
field  Is  decomposed  such  that  the  spatial  variations  of  the  field  depend 
on  wavelength.  As  written  In  Eq  (3-3)  the  aperture  field  Is  not  coherence 
separable.  If  for  Instance,  the  aperture  field  Is  due  to  a field  which 


had  propagated  some  distance  to  the  aperture  plane,  the  output  field 
consists  of  a sum  of  Individually  propagated  temporal  modes  (Ref  (2-42)), 
and  the  spatial  part  of  each  mode  is  given  by  Eq  (2-36) . To  simplify 


the  analysis  later  in  this  chapter.  It  will  be  assumed  at  a later  point 
that  the  Input  field  Is  coherence  separable,  and  as  such,  the  spatial 
variations  of  the  aperture  plane  are  the  same  for  all  temporal  modes. 
This  assumption  will  allow  the  calculation  of  the  phase  of  an  aperture 
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Figure  7 . Propagation  Geometry  for  X-Channel  Interferometer  Optics 

field  whose  spatial  content  is  Independent  of  wavelength.  If  U (r  ,t) 

a a 

were  coherence  separable,  Eq  (3-3)  would  be  written 

_ _ +j2iT^t 

for  (3-4) 

where  wavelength-independent  spatial  part  of  the  aperture 

field.  The  discussion  of  phase  measurement  for  a temporally  broadband 
field  which  is  not  coherence  separable  is  presented  in  Chapter  IV  where 
the  Interferometer  is  applied  to  the  situation  where  the  aperture  field 
has  a phase  which  has  been  distorted  by  atmospheric  turbulence. 

Using  the  free-space  propagation  model  for  temporally  broadband  fields 
expressed  in  Eq  (2-38) , the  n^^  temporal  mode  of  the  optical  field  Just 
prior  to  lens  L]  can  be  written 


( 

! 

( 
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Jk.Fi 


- T ^ T 
for  -2  <*^<2 


where 


U (r  ,t)  = u Cr  )w  e 
an  a'  an  a n 


+j2TrYt 


(3-5) 


(3-6) 


represents  the  n temporal  mode  of  the  aperture  field,  and  U«  (r«  ,t) 

■^In 

Is  the  n^^  mode  of  the  field  at  lens  plane  just  prior  to  thin  lens  Lj 

Apochromatic  Lens  Model.  Following  Goodman  (Ref  20:81),  a lens  per- 
forms a quadratic  phase  transformation  on  the  incident  field.  For  a 
white-light  incident  field  U£^(r^^,t)  an  apochromatic  lens  performs  a 
quadratic  phase  transformation  on  the  broadband  field  in  the  following  way: 


r T ^ T 

for  - Y b < Y 


(3-7) 


where  U»  (ro  ,t)  is  the  n^  temporal  mode  of  the  field  Just  to  the 
•<^ln 

right  of  lens  Lj  . The  n^**  field  mode  just  prior  to  the  grating  in 

plane  given  by  the  result  of  propagating 

U»  (r«  ,t)  a distance  Fi  to  the  focal  plane  of  the  lens: 

■^In  "^1 

jk  Fi 


««pCj2^|rg-r^,l^] 


(3-8) 


I 

IxTT 


I [-3  Pg-'g,] 


By  substituting  Eq  (3-5)  Into  (3-9)  and  solving  the  Integral  over 
assuming  the  lens  Is  of  Infinite  extent,  yields 


U 

gn 


.t) 


dr  U (r  ,t)exp[-j:=^  r • r ] 
an  a Fi  g a'' 


for 


T ^ T 


(3-10) 


Note  that  for  each  temporal  mode,  the  spatial  part  of  the  grating  field  Is 
an  exact  two-dimensional  Fourier  transform  of  the  Input  spatial  field: 


U 

gn 


Jk  2Fl 


" r - - - 

— dr  u (r  )exp[-j  2Tr(-r-^)*r  ] 
JA  Fi  I a an  a A Fi  a-* 

n ^ n * 


(3-11) 


jir?r  V"']  ■'xy 

n do 


(3-12) 


where  ^gn^’^g^  time- Independent  parts  of  the  aper- 

ture and  grating  fields,  respectively;  F {*)  denotes  a two- 

*a^a 


dimensional  Fourier  transform  with  respect  to  variables  x and  y ; 

a a 

and  the  transform  Is  evaluated  at  spatial  frequencies  f^  ■ 

and  fy  ■ • The  broadband  field  at  the  grating  Is,  therefore, 

the  superposition  of  Fourier  transforms  defined  for  each  temporal  mode. 

This  result  may  be  considered  a white-light  extension  to  the  Fourier 

transforming  properties  of  F-F  optical  systems  (Ref  20:86-87)  where  the 

Input  field  can  be  decomposed  as  In  Eq  (3-3) . If  the  lens  L;  were  not 

apochromatlc  as  defined,  the  transform  results  (3-11)  and  (3-12)  would  not 

hold  for  the  field  at  grating  plane  P . 

g 
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Rotating  Grating  Model.  Since  the  radial  grating  Is  assumed  to  be 

periodic  In  for  the  x-channel  of  the  Interferometer,  Its  transmission 

function  can  be  expanded  In  a Fourier  series.  The  n^^  mode  of  the  chopped 

field  just  beyond  the  rotating  grating  U'  (r  ,t)  can  then  be  written 

gn  g 


U^(r  ,t)  = U (r  ,t)  y G exp[j2Tr^(x  -v  .)] 
gn  g’  gn'  g’  ' i m d ' g gt'-* 


(3-13) 


th 


where  U (r  ,t)  Is  the  n temporal  mode  of  the  optical  field  Incident 
g g 

on  the  rotating  radial  grating  given  In  Eq  (3-12) 


„ 1 ^2™  .mw. 

sln(^) 


(3-14) 


th 


Is  the  Ronchl  (square-wave)  grating  coefficient  of  the  m diffracted 
order  for  all  Integers  m 0 , where  Gq  “ ^ 


2TrR 

N 


(3-15) 


Is  the  grating  period  of  an  N line  radial  grating  at  distance  R from 
the  center  of  the  disc,  and 


Re 


(3-16) 


Is  the  speed  of  the  grating  In  the  x -direction,  6 la  the  angular  speed 

g 

of  the  grating  In  radians  per  second. 

Note  that  the  phase  of  the  grating  coefficients  G given  In  Eq 

m 

(3-14)  Indicates  an  asymmetrical  positioning  of  the  grating  with  respect 

to  the  x -axis  at  t - 0 . Specifically,  a time  reference  may  be  defined 
g 

so  that  t ■ 0 whenever  the  position  of  the  grating  allows  Its  Fourier  coef- 
ficients to  be  given  by  Eq  (3-14) . This  time  reference  Is  necessary  for  the 
extraction  of  the  phase  measurement  from  the  Interferometer  output  signal. 


T 

i 

I 
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Detector  Plane  Field.  Using  the  broadband  propagation  result  for 

apochromatlc  lens  F-F  optical  systems  similar  to  the  one  expressed  In 

Eq  (3-10),  the  n^^  optical  field  mode  at  the  detector  plane  U,  (r,,t) 

dn  d 

due  to  the  broadband  field  U'(r  ,t)  Is 

g g 


, '•'g  'd-'g! 


r T ^ T 

for  - J < t < Y • 


(3-17) 


Combining  Eqs  (3-10),  (3-13),  and  (3-17)  and  carrying  out  the  two- 

dimensional  Integral  over  the  grating  plane  P yields  an  expression 

g 

for  the  n temporal  mode  of  the  detector  field  in  terms  of  the  aperture 
field: 


^dn^’^d’*^^  ■ I M exp[j2k^(Fi+F2)]exp[-j27Tmf^t] 


T T 

- TT  < t < ■;r 
2 2 


(3-18) 


where 


F2 


(3-19) 


is  the  magnification  of  the  two  lens,  double  F-F  optical  system,  f^ 

Is  the  fundamental  modulation  frequency  of  the  detector  field  given  pre- 
viously in  Eq  (3-2) 


(3-20) 
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Is  the  shear  distance  for  the  n temporal  mode,  and  sq  Is  the  shear 

distance  for  the  median  wavelength  Xg  , measured  In  the  detector  plane. 

The  expression  for  given  In  Eq  (3-6). 

The  detector  field  of  the  x-channel  for  each  temporal  mode  n and 

each  diffracted  order  m Is,  therefore,  a scaled  and  shifted  Image  of 

the  Input  field  mode  U (r  ,t)  modulated  at  frequency  mf  . Note  that 

dn  Si  s 

the  modulation  frequency  and  scaling  of  the  field  are  Independent  of  wave- 
length. The  location  (or  relative  shear)  of  each  diffracted  order,  however, 
depends  on  temporal  mode — and  therefore,  on  wavelength.  The  amount  of 
x-shear  for  each  diffracted  order  m Is,  from  Eq  (3-20) 


X Fo 
n ^ 

ms  • m — -z — 
n d 


(3-21) 


The  effect  of  this  wavelength  dependence  on  the  Interferometer's  phase 

measurement  for  aperture  fields  with  various  states  of  spatial  coherence 

Is  shown  In  the  last  two  sections  of  this  chapter. 

Following  the  propagation  model  for  broadband  fields  (2-40)  the 

broadband  detector  field  U (r.,t)  Is  written  as  the  sum  of  all  the 

a d* 

individually  propagated  modes: 


U^(r^,t)  - I U^^Cr^.t) 


(3-22) 


Analysis  of  Detector  Plane  Intensity  Pattern 

The  llght-lntenslty  on  the  detector  plane  l^(r^,t)  due  to  the 
broadband  field  Is  given  by 


n n ' f 


(3-23) 


(3-24) 
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where  n and  n'  range  over  all  significant  temporal  modes  of  the  broad- 
band process  w(t)  . Note  that  each  field  mode  given  in  Eq 

(3-18)  is  the  sum  of  all  the  field  orders  m diffracted  by  the  rotating 
grating.  The  interference  of  all  diffracted  temporal  modes  must  therefore 
be  written  as  a double  sum  over  m . The  detector  plane  Intensity  may 
also  be  written 


t,- 


”11  m2  G G exp[-j2TT(m-mOf„t] 


r_ (M(r .-ms  ) ,M(r  -m"s  ) ) 
a d n d n 


where 


r (r  ,r'  ) 
a a a 


^ ri  w 1 \ an  a an  a / 


(3-27) 


(3-28) 


(3-29) 


Is  the  spatial  correlation  of  the  broadband  aperture  field  U (r  ,t) 

a a 

evaluated  at  points  r and  r'  . 

3 SL 

From  Eq  (3-27)  it  can  be  seen  that  the  intensity  pattern  modulated 

in  time  at  frequency  f^  results  from  the  coherent  addition  of  two 

diffracted  orders  m and  m'  such  that  |m-m'|  = 1 . Using  the  Ronchi 

grating  coefficients  defined  in  Eq  (3-14) , the  Intensity  pattern 

modulated  at  f is  calculated: 
s 


" T Re{exp[-j(2irf^t-|)]r^(M(r^-a^),Mr^)} 


+ — Re{exp[+j(2iTfgt-y)]r^(M(rj+3^),Mr^)} 


T T 
for  - "2  ^ Y 


(3-30) 


where  >4^(r^,t)  is  the  intensity  pattern  at  the  detector  plane  at  fre- 
quency f^  for  x-shear,  and  the  double  sum  over  diffraction  orders  is 
calculated  such  that  Im-m'l  •*  1 for  all  terms  in  the  summation.  The 
intensity  pattern  .(.^(r^.t)  at  each  point  (Xjj»y^)  the  detector 
plane  is  defined  to  be  the  x-channel  output  signal  of  the  AC  interfero- 


meter. 
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Each  detected  signal  mode  n Is  composed  of  two  similar  terms  modu- 


lated at  frequency  f . The  first  term  is  the  result  of  the  coherent 

s 

addition  of  the  m • +1  and  m = 0 field  orders  diffracted  by  the  rotating 
grating.  This  term  will  be  called  the  "positive-shear"  term  because  it 
contains  the  correlation  of  the  +1  diffracted  order,  sheared  by  distance 
+8^  in  the  detector  plane,  with  the  undiffracted  field  order.  The  second 
term  results  from  the  interference  of  the  m ••  -1  and  m = 0 diffracted 
orders  and  is  called  the  "negative-shear"  term. 

Since  the  spatial  part  of  the  aperture  field  mode  is  Independent  of 
time,  no  additional  frequency  modulation  is  introduced  into  the  detector 
signal. 


( 


Phase  Measurement  for  Coherence  Separable  Aperture  Fields 

In  this  section  the  Interferometer  phase  measurement  is  derived  for 
aperture  fields  which  may  be  assumed  to  be  coherence  separable.  Such 
fields  which  can  be  characterized  by  an  expansion  over  its  temporal  modes 
may  be  written  as  in  Eq  (3-A) . Following  Eq  (3-28) , the  spatial  correla- 
tion of  the  coherence  separable  aperture  field  may  be  written 


where 
field 
to  be 


(3-31) 

(3-32) 


n 

R (r  .r")  is  the  correlation  of  the  spatial  part  of  the  aperture 

d fi  d 

defined  in  Eq  (2-5).  Allowing  the  spatial  correlation  R (r  ,r') 

& A Q 

written  in  terms  of  its  complex  degree  of  coherence  (Ref  Eq  (2-53)) 


R (r  ,r'  ) - .r'  )u  (r.)u*(r'  ) 

aaa  aaa  aaaa 


(3-33) 
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where  0 (r  ,r')  is  the  complex  degree  of  coherence  of  the  aperture  field, 

A & & 

the  detector  signal  for  x-shear  may  be  written 


2 

-tjCrd.t)  « ^ Re{exp[-j  (2Trf^t  - -j)]  u^(Mr^) 

n 

+ ^ Re{exp[+j(2iTf^t 

n 


for  (3-34) 

where  all  the  mode-dependent  terms  have  been  grouped  together. 

Let 


and 


( 


V 


S (^)fl  (M(r.-s  ) ,Mr  ) 
w T a and 


. u (M(r  .-s  )) 
a d n 


(3-35) 


ujM(r^+so))  I 


. u (M(r.+8  )) 
a d n 
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(3-36) 
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Furthermore,  let  all  apatlal  fields  be  written  In  polar  form: 


Ua(Mr^)  - A(Mr^)exp[j<t.(Mr^)  ] 


(3-37^ 


u^(M(r^-so))  = A(M(r^-so))exp[j(fi(M(r^-So))] 

and 


(3-38) 


Ua(M(rd+8o))  “ A(M(rj+so))exp[j(J.(M(r^+8o))] 


(3-39) 


( 


In  the  next  section  it  will  be  shown  that  for  white-light  coherence 
separable  aperture  field  envelopes,  the  lateral  shear  AC  interferometer 
causes  the  m = +1  and  m * -1  diffracted  spatial  field  orders  to  be 
spatially  low-pass  filtered.  The  "hats"-  over  the  amplitude  and  phase  of 
the  diffracted  orders  indicate  that  they  represent  the  amplitude  and  phase 
of  the  spatially  filtered  field.  The  nature  of  this  spatial  filtering  is 
generally  discussed  for  coherence  separable  broadband  fields  in  the  next 
section,  and  discussed  specifically  for  broadband,  spatially  Incoherent 
source  fields  in  Chapter  IV.  Note  that  the  spatially  filtered,  diffracted 
field  orders  are  defined  vlth  respect  to  the  shear  distance  sq  for  the 
median  wavelength  of  white-light.  Whereas  the  performance  of  the  inter- 
ferometer is  analyzed  with  respect  to  sq  , it  may  be  used  as  a design 
parameter  for  a physically  realizable  Interferometer. 

Note  that  the  m > 0 , or  undiffracted  order.  Is  not  filtered.  The 
spatial  filtering  of  the  diffracted  orders  can  be  thought  of  as  being 
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I caused  by  differential  path  length  effects  associated  with  the  diffrac- 

tion of  broadband  light  by  the  Interferometer's  rotating  grating.  If  the 
field  Is  not  diffracted,  there  are  no  differential  paths  and  no  spatial 
filtering  Is  Introduced. 

Using  the  expressions  given  above  for  the  zero  (3-37),  m » +1, 
(3-38),  and  m ■ -1  (3-39)  diffracted  field  orders,  Eq  (3-34)  may  be 

rewritten 

'td(7^*t)  - ^ A(M7^)A(M(7^-Io))Re{exp[-j(27rf^t-|  )] 
exp  [ j (KM(  r ^-s  0 ) ) - j ()•  (Mr^)  ] } 

u2  ^ ^ ^ 

+ — A(Mr^)A(M(r^+so)Re{exp[+j(27rf^t -y  )] 

( - - - ' 

exp  [ j (j)  (M(r^+s  o ) ) (Mr^)  ] ) 


( 


T T 

for  ” 7 ^ ^ "2  • 


(3-40) 


If  the  amplitudes  of  the  spatially  filtered  orders  are  assumed  to  be 
equal  over  the  distance  of  the  shear  for  each  point  In  the  detector  plane 


A(Mrj)  = A(M(r^-So)) 


“ A(M(r,+so)) 


(3-41) 

(3-42) 


where  A(x^,y^)  Is  the  amplitude  of  both  sheared  field  orders,  the  posi- 
tive and  negative-shear  terms  may  be  combined  and  Eq  (3-40)  may  be 
simplified: 
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'Cd(rd,t)  - — jj-  A(Mrd)A(Mrd) 


cos  [-j  )-'!'] 


8in[2TTf^t+-|-  ] 


- T ^ T 

for  “ "2  t <-2 


(3-A3) 


where 


- (MTd) 


(3-4A) 


♦+  - -KMCrd-sJ) 


(3-45) 


- (j.(M(rd+8o)) 


(3-46) 


Eq  (3-43)  represents  the  interferometer  output  signal  for  x-shear  in  terms 
of  the  amplitudes  and  phases  of  the  sheared  (and  filtered)  aperture  field 
orders.  The  aperture  phase  measurement  made  by  the  Interferometer  output 
signal  Is 


Ai^  * ^ - 


(3-47) 


It  is  shown  In  later  sections  how  the  measured  phase  A4i  relates  to  the 
spatial  part  of  the  broadband  aperture  field  for  special  cases  of  aperture 
field  coherence. 

If  the  amplitudes  of  the  spatially  filtered  orders  are  not  equal  over 
the  distance  of  shear  as  assumed  In  Eqs  (3-41)  and  (3-42),  the  resulting 
detector  signal  is  rewritten 
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^d(rd,t) 


m2  — 

- ^ A(Mr .) 

TT  a 

{A(M(rd-So)  ) sln[2Trf ^t-(()(M(rd-So) )+ij>(Mrd)  ] 

+ A(M(r,+So))  sin[2Trf  t 
d “ s 

+ ♦(M(7d-*^o))-'KM7d)]} 


(3-48) 

Eq  (3-48)  can  be  Interpreted  as  a sum  of  two  phasors  with  phases 

>p  - -<KM(rd-so))  + 4>(M7d)  (3-49) 

due  to  the  positive-shear  term,  and 

- +i(M(7d+7o))  - 'KM7d)  (3-50) 

due  to  the  negative  shear  term  as  indicated  in  Figure  8.  The  resulting 

A A 

amplitude  A^  and  phase  of  the  detector  signal  is  the  sum  of  the 

positive-  and  negative-shear  phasor  terms  as  Indicated  in  Figure  8. 
Henceforth,  it  will  be  assumed  that  the  amplitudes  of  the  spatially  fil- 
tered diffracted  fields  are  nearly  equal  so  that  Eq  (3-43)  adequately 
represents  the  interferometer  output  signal  for  x-shear.  It  is  shown  in 
Chapter  IV  that  for  the  case  of  distant  sources,  this  assumption  is, 
indeed,  valid. 
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If  D Cr  ifp  Is  homogeneous  across  the  aperture  of  the  Interferometer, 

cl  & a 

the  complex  degree  of  coherence  of  the  aperture  field  may  be  written 


a 

a 


(3-52) 


Assuming  p homogeneous  allows  (3-51)  to  be  written 

Sl 


u (M(r  -so))  = ? I S„(^)P^  (Ms„)u  (M(r,-8„))  (3-53) 


T ^ ■'w'T"‘a '““n''*a'“'‘d  “n' 
n 


Note  that  p Is  evaluated  for  a separation  of  Ms  In  the  x-dlrectlon 
a n 

and  no  separation  In  the  y-dlrectlon.  Taking  the  spatial  Fourier  transform 
of  Eq  (3-53)  with  respect  to  the  aperture  plane  coordinates  yields 

^x  y 
a a 


1 n _ -j2Trf  Msor^(n) 

Y I )Wa(M8^)  V^(f^,fy)e 


where 


V (f  ,f  ) •«  F {u  (x  ,y  )} 
a X y X y a a ■^a 
a'^a 


(3-54) 

(3-55) 


Using  the  binomial  expansion  for  C~^(n)  and  approximating  the  Infi- 
nite sum  by  the  first  two  terms 


r^n)  . (1-^)-^  = 1^^  (3-56) 

which  Is  a good  approximation  for  « 1 , allows  (3-54)  to  be 

lOT 

written  In  the  following  form: 
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y ^u^<M(r^-so))}=  )exp[-j2Trf^Mso] 


exp[-j 2tt^x^^°  (^)] 

fo 


(3-57) 


The  approximation  In  (3-56)  assumes  that  although  the  power  spectrum 
S^(f)  Is  broad.  Its  effective  bandwidth  B Is  still  sufficiently  small 
that  the  higher  order  terms  of  C \n)  In  the  exponential  of  Eq  (3-5A) 
can  be  neglected.  The  maximum  Index  n In  the  sum  (3-57)  Is  limited  by 
the  effective  bandwidth  B so  that  Max(n)  = + BT  . Assuming  that  the 
power  spectrum  Is  still  broadband  enough  that  the  sum  over  n may  be 
approximated  by  an  Integral  over  v 


^x  y ^'^a(M(rj-so))  } ' V^(f^,f  )exp[-j2Ttf^Mso] 
a'^a  ^ 


jdv  q^(M8o(l+-^))  S^(v) 


f Mso 

exp[-j2Tr(  -2 ) v] 

to 


(3-58) 


To  first  order  then,  Eq  (3-58)  expresses  the  spatial  filtering  effect  of 
the  AC  shearing  Interferometer  on  the  m ■ +1  diffracted  field  order, 
when  the  source  field  can  be  assumed  to  be  coherence  separable.  Although 
Eq  (3-58)  Indicates  a doubly-lnfinite  range  of  Integration,  the  range  of 


53 


V for  which  the  integral  has  significant  value  is  limited  to  the  spectral 
bandwidth  of  S^(v)  — which  is  nominally  -B  < v < B . Define 

* 

= dvu^(Mso(l -■^))  S^(v)exp[-j2TTCv]  (3-59) 

where  C = (3-60) 

fo 

to  be  the  spatial  filter  introduced  by  the  interferometer's  diffraction 
grating  on  the  m = +1  field  order.  Note  that  there  is  no  spatial 
filtering  for  spatial  frequencies  associated  with  the  y-dlrection  since 
no  y-shear  was  introduced  in  the  x-channel.  The  corresponding  filter 
for  the  m = +1  order  in  the  y-channel,  therefore,  does  not  affect 
the  x-direction  spatial  frequencies. 

The  same  analysis,  when  applied  to  the  m = -1  diffracted  order, 
yields 

a^a 

- V^(f^,fy)exp[+j2Trf^Mso]H_(0  (3-61) 

where 

H (C)  - [ dv  0 (M8o(1+-^))  S (v)exp[+j2TTCv] 

J . *0  w 

is  the  spatial  filtering  Introduced  into  the  m ■ -1  diffracted  field 
order . 


( 


1 


The  form  of  spatial  filtering  functions  and  H is  seen  to 

depend  on  three  effects:  (1)  the  complex  degree  of  coherence  of  the 


aperture  field  u (r  -r')  , (2)  the  power  spectrum  of  the  wide-hand 

Si  3i  ^ 

temporal  process  S (f)  , and  (3)  the  shear  distance  sn  for  the  median 

V 

light  wavelength,  defined  by 


( 


So 


Mi 

d 


(3-63) 


The  shear  distance  is  a design  parameter  which  may  be  chosen  to  fit 

a specific  set  of  operating  conditions.  Since  the  coherence  properties 

of  broadband  light  are  relatively  insensitive  to  the  precise  functional 

form  of  the  temporal  process  power  spectrum,  the  ultimate  forms  of  the 

spatial  filtering  functions  are  also  relatively  insensitive  to  the  precise 

form  of  the  broadband  power  spectrum  S (f)  . The  dependence  of  the  fil- 

w 

terlng  effect  on  the  specific  coherence  properties  of  the  aperture  field 
is,  however,  significant.  In  the  next  two  sections,  the  effect  of 
aperture  field  coherence,  given  a particular  temporal  power  spectrum, 
are  considered  for  special  cases  of  the  field's  degree  of  coherence 
{]  (r  ,r')  . The  filtering  effects  of  a specific  form  of  p (r  ,r')  are 

d & S Add 

considered  In  Chapter  VI  with  the  application  of  the  shearing  inter- 
ferometer as  a wavefront  sensor  in  a phase-compensated  Imaging  system. 


Filtering  Effects  for  Spatially  Coherent  Aperture  Fields 

For  an  aperture  field  which  is  spatially  coherent,  the  complex  degree 
of  coherence  is  equal  to  unity  for  all  points  of  interest  in  the  aperture 
plane  (Ref  Eq  (2-34)).  Following  Eq  (3-59)  the  spatial  filter  may 
then  be  written 


\ 
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(3-64) 


H^(4)  ■ |dv  S^(v)exp[-j2TTCv] 

" F^{S^(v)}  (3-65) 

f Mso 

“ R ° -? > (3-66) 

W wig 

where  F^{ • } denotes  the  Fourier  transform  with  respect  to  variable  v , 
and  the  transform  is  evaluated  at  , and  R (•)  is  the  corre- 

—z W' 

latlon  of  the  broadband  time  process  w(t)  defined  in  Eq  (2-9) . Note 
that  since  the  correlation  functions  considered  here  are  stationary 
and  even 

F^{Sw(v)}  = f;;^  {S^(v)}  = R^(0  (3-67) 

Applying  the  same  analysis  to  Eq  (3-63)  shows  that  the  spatial 
filtering  Introduced  into  the  -1  diffracted  field  order  has  the  same 
form: 


H (C)  = R (^)  (3-68) 

~ w 

- H_^(c)  - H(c)  (3-69) 

As  discussed  in  Chapter  II  (Ref  Eq  (2-73)),  the  correlation  function 
R^(t)  will  drop  off  with  increasing  argument  so  as  to  limit  the  high 
spatial  frequency  content  of  V (f  ,f  ) in  Eq  (3-58).  Therefore,  for 

A X ^ 

tdilte-light,  spatially  coherent  aperture  fields,  the  AC  shearing  inter- 
ferometer causes  the  sheared  orders  of  the  broadband  detector  field  to 
be  spatially  low-pass  filtered.  The  form  of  this  filtering  has  been 
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A 


shown  to  be,  to  first  order,  the  correlation  function  of  the  broadband 


temporal  process  w(t)  evaluated  at  t = x 


f-Mso 


As  an  example  calculation,  consider  the  bandllmlted  power  spectrum 
Sj^(f)  and  its  correlation  function  R^(t)  = F^^{S^(f)}  . The  correla- 

tion function  resembles  a sin(x)/x  function  with  the  first  zero-crossing 
at 


fxMso 


2B  "^b 


(3-70) 


where  B is  the  bandwidth  of  S^(f)  . is  the  argument  of  Rjj(t) 
for  which  ^(t)  =0  , and  the  argument  of  the  spatial  filtering  func- 
tion H(c)  has  been  equated  with  that  of  the  correlation  function 
Define  the  effective  spatial  bandwidth  of  the  filtering  function 
H(  ^x^^°  ) for  spatially  coherent  aperture  fields  to  be 

fo 


B, 


H 2MsnB  Msn  ' Af 


Af 


(3-71) 


For  white-light  applications,  =>  0.275  and  the  effective  spatial 
bandwidth  of  the  spatial  low-pass  filter  H(c)  is  on  the  order  of 


3.636 

Msq 


(3-72) 


or,  in  other  words,  the  spatial  frequency  content  of  the  sheared  broad- 
band fields  u^(M(r^-So))  and  u^(M(rj+so))  is  limited  to  3.6  times 
the  reciprocal  shear  distance  Msp  measured  in  the  aperture  plane. 
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since  the  spatial  field  u (r  ) Is  complex,  HC^)  filters  both 

Si  d 

real  and  Imaginary  quadratures  of  the  sheared  fields.  The  phase  measure- 
ment derived  from  Eq  (3-43),  however.  Is  expressed  In  terms  of  the 
magnitude  A and  phase  <j)  of  the  filtered  spatial  part  of  the  aperture 
field  envelope.  The  specific  form  of  the  distortion  Introduced  Into 
the  measured  aperture  phase  (p  by  the  filtering  process  cannot.  In 
general,  be  determined  for  arbitrary  forms  of  ()i(r  ) . At  this  point. 

It  will  suffice  to  say  that  the  phase  of  a white-light  field  measured  by 
the  AC  shearing  Interferometer  Is  the  phase  of  a filtered  aperture  fleld- 
the  filtering  being  Introduced  by  diffraction  of  the  white-light  field 
by  the  Interferometer's  grating. 

Monochromatic  Field  Phase  Measurement.  For  spatially  coherent  aper 

ture  fields  which  are  monochromatic,  the  correlation  function  R (t) 

w 

Is  flat  so  that  H(^)  Is  constant  for  all  spatial  frequencies.  For  the 
monochromatic  case,  then,  no  spatial  filtering  Is  Introduced  for  any 
shear  distance  Msg  • With  no  spatial  filtering,  the  resulting  Inter- 
ferometer output  signal  for  x-shear  Is  given  by 

^^(7^,t)  » M2[A(r^)]2  cos[-| 

8ln[2irf^t -y  (())_-<|)^)  ] (3-73) 

where  A(r  ) Is  the  amplitude  of  the  aperture  field  at  point  r ■ Mr. 

a ad 

♦ - <K7^)  (3-74) 

Is  the  phase  of  the  spatial  part  of  the  aperture  field 

■ ♦(r^-Msg)  (3-75) 
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and 


♦ - <C(r  -mso)  (3-76) 

— a 

Note  that  for  aperture  phase  functions  (t)(r  ) which  change  at 

s 

most  linearly  at  each  point  in  the  aperture,  the  detector  signal  may 
be  simplified  further: 

'i-d(r^.t)  - ^ M2[A(7^)]2sin(2Trf^t-|A4,)  (3-77) 

where 


A()i  ■=  (J((r  +Mso)  - i(i(r  -Msq) 

a SL 


(3-78) 


is  the  linear  change  in  phase  over  distance  Ax^  « 2Mso  measured  in  the 
aperture  plane  for  a monochromatic,  spatially  coherent,  aperture  field. 
Referlng  to  Figure  9,  the  local  wavefront  slope  at  a point 


ao 


(x^^,y^^)  in  the  aperture  plane  is  approximately  given  by 

*<Xa/Mso,ya^)  - »(x^„-Mso  ,y^^) 


T — 4>(x  ,y  ) 

3x  a”ao 


X *x 
a ao 


2Msf 


(3-79) 


A(t> 
2M8  0 


(3-80) 


The  validity  of  the  approximation  depends  on  how  much  the  aperture  phase 
actually  fluctuates  within  a distance  Ax  ■ 2Mso  . A detailed  study  of 

3i 

aperture  phase  fluctuations  based  on  the  phase  statistics  of  turbulent 
atmosphere  can  be  found  in  Ref  24. 
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♦(x  .y„  ) 


(|i(x  +Mso,y^  ) 


♦(x^^-Mso.y^^)  4^- 


X -Msn 

ao 


X +M8o 

ao 


Figure  9.  Aperture  Phase  (|)(x  ,y  ) for  Fixed  y Showing  Calculation 
of  Local  Wavefront  Slope  a ao  a 


Therefore,  for  monochromatic  light,  the  slope  of  the  aperture  field 
wavefront  Is  given  by  Eq  (3-79).  Implicit  assumptions  are  that  the  aper- 
ture field  amplitude  Is  constant,  and  the  phase  of  the  aperture  field 
varies  at  most  linearly  over  a distance  2Mso  aperture  plane. 

Phase  Measurement  for  Small  Shear.  For  spatially  coherent,  white- 
light  aperture  fields  where  the  shear  Mso  Is  small  enough  that 

R ( ^x"°°  ) encloses  all  the  spatial  frequency  components  of  “ (r^) 
fo 


“s  ‘ ^ 


(3-81) 


where  B Is  the  bandllmltlng  spatial  frequency  of  V (f  ,f  ) , the 
A A X y 

spatial  filtering  effect  of  the  shearing  Interferometer  Is  negligible. 

A A 

For  this  case  the  filtered  phases  and  defined  In  (3-45)  and 
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(3-46)  very  closely  approximate  the  actual  phases  and  <J)  of  the 

sheared  Images  of  the  aperture  field  envelope.  The  detector  signal, 
then.  Is  expressed  by  Eq  (3-73)  and  the  wavefront  slope  for  x-shear  Is 
expressed  by  (3-80) . 

Phase  Measurement  for  Large  Shear.  For  spatially  coherent,  white- 

f Ms 

light  aperture  fields  were  the  shear  Msn  Is  large  so  that  R ( x^°  ) 


is  much  narrower  than  the  spatial  frequency  content  of  u (r  ) 

a a 

“a  ''  'h  " ^ 


(3-82) 


the  phase  measurement  A({i  Is  lost  completely.  Using  Eq  (3-43)  and 
approximating  R ( ) by  an  Impulse  at  f *0  In  Eq  (3-66)  yields 

W ^ X 

to 


where 


^d(rd*t)  “ M^AoA(r^)cos(({i-i{io)sin(2Trf^t) 

n 

B » — ^ (■^) 

H Msq  '■  Af 

Ao  = ijdx^  u^(x^,y^)| 


(3-83) 


(3-84) 


(3-85) 


is  the  amplitude  of  the  aperture  field  averaged  over  x , and 


- arg[  |dxau^(x^,y^)] 


(3-86) 


is  the  phase  of_  the  aperture  field  averaged  over  x 


Filtering  Effects  for  Aperture  Fields  with  Arbitrary  Spatial  Coherence 
For  broadband  aperture  field  envelopes  with  a complex  degree  of 
coherence  , the  first-order  spatial  filter  for  the  m ■ +1  diffracted 


field  order  Is  given  by  Eq  (3-59) : 
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f Mso  „ 

H.  ( -7 )-  dvO  (^p(fo4v))S  (v) 

+ fg  J a fg  w 

f Msg 

exp[-j27r(-|^)  v]  (3-87) 

Note  that  although  p has  arbitrary  form.  It  is  assumed  to  be  stationary 

d 

for  all  separations  r -r'  In  the  Interferometer  aperture  plane. 

d d 

The  spatial  filtering  function  can  be  Interpreted  as  a Fourier 

transform  of  a "windowed"  version  of  the  temporal  power  spectrum: 


f Msg 

H - F^{S  (v)}  (3-88) 

+ fg  V w 

where 

S„(v)  - q(^(fg+v))S Jv)  (3-89) 

W iQ  W 

Is  the  windowed  version  of  the  power  spectrum,  and  the  transform  Is 

evaluated  at  f^Msg/fg  . The  windowing  effect  can  be  better  visualized 

with  the  help  of  Figure  10,  where  the  windowing  Is  accomplished  by  a 

shifted  copy  of  the  complex  degree  of  coherence  p as  prescribed  by 

a 

Eq  (3-87) . The  degree  of  coherence  sketched  In  Figure  10  Is  equal 

to  the  modulus  of  the  complex  degree  of  coherence  p (Ref  14:510). 

a 

The  spatial  filtering  function  , then.  Is  the  Fourier  transform 
of  the  windowed,  or  filtered,  power  spectinim  sketched  In  Figure  10(b). 

Similarly  following  Eq  (3-62) , the  spatial  filter  for  the  m ■ -1 
diffracted  order  Is  the  Inverse  Fourier  transform  of  the  windowed  power 

A 

spectrum  S^(v)  defined  by  Eq  (3-89) 

f Msg 

H^(-^)  - f‘^{S^(v)}  (3-90) 
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Figure  10.  Windowing  Effect  of  Degree  of  Coherence  Pg  on  the 

Temporal  Power  Spectrum  : (a)  Power  Spectrum  and  Degree  of 

Coherence,  (b)  Windowed  Temporal  Power  Spectrum  S (v) 

w 

It  can  be  seen  from  Eqs  (3-88)  and  (3-90)  that  the  specific  form  of 
the  filtering  processes  and  H are  highly  dependent  on  the  complex 

degree  of  coherence  of  the  aperture  field.  Whereas  the  spatial  filtering 
arises  from  the  broadband  temporal  characteristics  of  the  field,  the 
functional  form  of  this  filtering  depends  on  the  spatial  coherence  of  the 
field  as  well. 

In  the  next  chapter  the  lateral  shear  AC  interferometer  is  applied 
to  the  measurement  of  an  aperture  phase  due  to  a broadband,  incoherent 
source.  It  is  found  that  for  specific  source  radiance  distributions  the 


the  complex  degree  of  coherence  of  the  aperture  field  can  be  calculated 
explicitly,  and  the  particular  form  of  the  spatial  filtering  functions, 
for  aperture  fields  which  are  not  assumed  to  be  coherence  separable  can 
be  calculated. 


IV.  WAVEFRONT  MEASUREMENT  FOR  SPATIALLY  INCOHERENT, 

TEMPORALLY  BROADBAND  SOURCES 

The  lateral  shear  AC  Interferometer  described  In  Chapter  III  has  been 
used  as  a wavefront  sensor  In  real-time  atmospheric  compensation  systems 
(Refs  7,  8,  9,  and  10).  In  this  application,  the  Interferometer  Is  required 
to  measure  the  phase  of  the  optical  field  envelope  In  Its  Input  aperture 
due  to  an  extended,  white-light  source.  This  chapter  describes  the  closed- 
loop  operation  of  the  lateral  shear  AC  interferometer  as  a wavefront  sensor 
for  broadband,  spatially  Incoherent,  extended  sources  based  on  the  free- 
space  propagation  model  Introduced  In  Chapter  III.  The  residual  effects 
of  atmospheric  turbulence  are  modeled  as  a multiplicative  phase  factor  Intro- 
duced at  the  aperture  plane  of  the  wavefront  sensor  for  each  temporal  mode 
of  the  field  expansion  for  a frozen  state  of  turbulence. 
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Figure  11.  Free-Space  Propagation  Geometry  for  Extended,  White-Light 
Source 

where  u^(r^)  la  the  spatial  part  of  the  source  field  envelope,  the  sum 
over  the  temporal  modes  converges  In  mean-square. 

If  the  source  Is  sufficiently  diffuse,  the  light  reflected  by  It  can 
be  assumed  to  be  spatially  Incoherent;  the  spatial  correlation  of  the  source 
field  may  be  written  (Ref  24:126  and  Eq  (2-55)) 


Aperture  Field.  Using  the  free-space  propagation  model  for  broadband 
fields  developed  In  Chapter  II  (Ref  Eq  (2-4A)),  the  n^*'  temporal  mode  of 
the  optical  field  U^Cr^.t)  just  prior  to  the  aperture  plane  is  given 

by 


. U (r,t)*u  (r)we 
pn  a’  an  a n 


+j2lT^t 


T T 

for 


where 


jk  Z 


“an^^a^  “ jA  Z 


" r _ 


(4-3) 


(4-4) 


Note  that  since  the  field  Up(r^,t)  is  the  output  field  due  to  a spatially 

Incoherent,  coherence  separable,  broadband  source  field  U (r  ,t)  , the 

8 8 

spatial  correlation  of  the  field  Is  given  by  Eq  (2-57) ? 


r (r  ,r'  ) 
p a’  a 


(r  .t)U*(r' 
pa  P n 


^ XgZ  ^ I T T 


j(l> 

n \ - n 
w'  T 


(4-5) 


(4-6) 


where 


(4-7) 
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Turbulence-Induced  Phase  Model.  The  effects  of  the  atmosphere 
seriously  degrade  the  imaging  capability  of  optical  systems  in  many 


applications.  The  resolution  which  is  attainable  for  imaging  through 
the  earth's  atmosphere  is  limited  by  warping  of  the  Isophase  surfaces 
(phase  distortion)  and  intensity  variations  across  the  wavefront  caused 
by  random  fluctuations  in  the  refractive  Index  of  air.  For  broadband 
optical  fields,  the  phase  distortion  which  is  introduced  can  be  thought 
of  as  a differential  path  length  distortion  introduced  for  all  wavelengths, 
so  that  the  phase  distortion  for  each  temporal  mode  measured  in  the  aper- 
ture plane  of  the  interferometer 


I 


( 


» 

) 

» 


n 

where  AZ(r  ) is  the  differential  path  length  distortion  introduced  by 
a 

a frozen  atmospheric  state  for  all  temporal  modes  measured  at  a point  r^ 
in  the  aperture  plane. 

To  the  extent  that  the  phase  fluctuations  correlated 

over  the  entire  wavelength  range  AX  of  the  broadband  source,  the  phase 
distortion  for  each  temporal  mode  is  approximately  the  same,  and 


for  -BT  < n < BT  (4-9) 

irfiere  ^ (r  ) is  the  turbulence-induced  phase  distortion  for  X * Xq, 
Tq  a 

which  is  the  same  for  all  wavelengths  of  light.  For  arbitrary  spectral 
bandwldths,  however,  the  phase  distortion  may  not  be  assumed  to  be 
correlated  for  all  modes  of  the  field  expansion  (Ref  28:468). 
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I ^ 


Suppose  that  the  mode  phase  distortion  is  written  as 


(4-10) 


where  <1'™  (r  ) is  the  fluctuation  of  the  turbulence- induced  phase  for 

m,  Xu  d 

each  temporal  mode  about  a mean  phase,  which  is  arbitrarily  chosen  to 

be  that  for  Xn-  In  general,  ({i_  (r  ) will  have  values  which  differ 

_Tn  a 

greatly  from  'll-,  (r  ) . If  gross  differential  path  length  differences 

IQ  3. 

and  wavefront  tilt  have  been  removed, as  in  the  closed-loop  operation  of 

a real-time  atmospheric  correction  device  (Refs  7 and  9),  the  variation  of 

the  phase  perturbation  about  (|i_  (r  ) for  each  mode  will  be  about  zero- 

Tq  a 

mean  over  the  range  of  the  modal  expansion: 


for  -BT  < n < BT 


(4-11) 


where  the  expected  value  may  be  taken  for  each  turbulent  state  of  the 
atmosphere.  Furthermore,  if  attention  is  restricted  to  points  in  the 
aperture  plane  such  that 

Max[|r^-r'  |]  < Po  (4-12) 

where  pg  is  the  spherical-wave  coherence  length  for  the  space-to-earth 
propagation  path  (Refs  2:1376  and  26:1678),  the  correlation  of  the  phase 
distortion  for  each  monochromatic  mode  of  the  field  expansion  may  be 
written  (Refs  26:1683;  27:731;  and  28:464-468) 

This  result  Is  useful  In  the  evaluation  of  the  Interferometer  phase 
measurement  for  restricted  operating  conditions. 
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For  real-world  turbulence  conditions,  the  random  phase  (i  (r  ) 
fluctuates  with  time  causing  scintillation.  In  order  for  the  phase  com- 
pensation system  to  operate  properly,  the  response  time  of  the  Imaging 
system,  and  so  the  wavefront  sensor,  must  be  short  enough  to  respond  to 
significant  temporal  fluctuations  In  the  phase  distortion.  This  require- 
ment limits  the  measurement  Interval  of  the  Interferometer — and  also  the 
characterization  Interval  of  the  broadband  process  T — tc  less  than 
about  10  milliseconds  (Refs  26:1680  and  29:392).  Note  that  this  require- 
ment does  not  prohibit  the  use  of  the  modal  propagation  model  since  any 
practical  measurement  Interval  Is  still  much  greater  than  the  reciprocal 

bandwidth  of  the  temporal  fluctuations  of  the  white-light  envelope 
-13 

(about  6.7  X 10  seconds). 

Detector  Plane  Output  Signal.  Modeling  the  effects  of  atmospheric 
turbulence  as  a unit-modulus  phase  screen  for  each  temporal  mode,  the 
n temporal  mode  of  the  field  In  the  Interferometer's  aperture  plane 
can  be  written 


(4-14) 


where  ” mode  of  the  aperture  field  due  to  an 

Incoherent,  broadband  source  given  previously  In  Eq  (4-3). 

Writing  the  broadband  aperture  field  as  the  superposition  of  all  temporal 


modes  yields 


U (r  .t)  - I u^(r  )w  e T ^ Tn  a 

urn  ^ An  A n 

n 


for 


(4-15) 


Assuming  that  the  random  phase  fluctuations  of  the  turbulent  atmo- 
sphere are  Independent  of  the  source  statistics,  the  spatial  correlation 


of  the  aperture  field  can  be  calculated  as  follows: 
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( 


(4-16) 


'■.‘'.•'I  > ■ •']> 

■ I ^ 

for  -•|<t<|  (4-17) 

>*«'  <$an'V"L<';)  ■<17>4vf> 

f _ _ k 

Note  that  for  closed-loop  operation  of  an  atmospheric  correction  system, 
the  statistics  of  the  phase  fluctuations  and  the  modal  expansion  coeffi- 
cients not  be  completely  uncorrelated.  It  Is  assumed  that,  for 

normal  operation  of  the  Interferometer,  the  correlation  of  these  two 
quantities  Is  small  enough  that  the  correlation  of  the  aperture  field  can 
be  written  as  In  Eq  (4-17) . 

Without  loss  of  generality,  the  turbulence-induced  phase  perturba- 
tion for  each  temporal  mode  may  be  represented  as  In  Eq  (4-10) , and 
Eq  (4-17)  may  be  rewritten: 

[ ^ (r  )u  (r'  N 
“ \ an  a an  a / 

^xp[j^Tn(ra)-J^Tn<*^;  )]N 

for  • (4-19) 

( 
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For  the  case  where  the  phase  fluctuations  correlated  over 

the  temporal  bandwidth  of  the  source,  (p-  (r  ) Is  Identically  zero  and 

_Tn  a 

the  last  expectation  Involving  "^Tn^^a^  will  yield  unity. 

For  a real-time  atmospheric  correction  system  operating  In  closed-loop, 

(J^(r^)  represents  the  residual  phase  error  associated  with  tracking  the 

average  phase  <|>„  (r  ) of  the  aperture  field, 
i 0 ^ 

As  an  example  calculation,  assume  the  closed- loop  system  Is  tracking 
well  so  that  zero-mean  (Ref  Eq  (4-11)).  Furthermore,  assume 

^^(r^)  Is  normally  distributed  so  that  the  last  expectation  In  Eq  (4-19) 
can  be  written  In  terms  of  the  characteristic  function  of  a Gaussian  ran- 


dom variable  (Ref  30:159-160): 


where 


^xp[j^Tn^’^a>-J!Tn^<  )]^  = 

W^a' 


(4-20) 


(4-21) 


If  the  statistics  of  (r  ) are  stationary  In  P , the  variance  of 

_Tn  a a 

the  residual  phase  fluctuations  may  be  further  simplified: 


(4-22) 


wbere 


^Tn^’^a^tln^’^a  V 


^1—  \-in  a -in 


(4-23) 


ru  ilisiTs.  If  the  phase  measurement  Is  restricted  to  regions  In  the 
«M»ae«  aa«lafTla«  Eq  (4-12),  nearly  one  for  all  points 

•«  ami  CIm  eapactatloo  trrlLten  In  Fq  (4-20)  Is  nearly  unity. 
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The  above  discussion  has  Intended  to  demonstrate  that  for  the 


* 


restricted  cases  discussed  above,  namely  (1)  the  phase  fluctuations 
♦,j^(r^)  are  correlated  over  the  entire  wavelength  range  of  the  source 
emission  and  (2)  the  residual  phase  error  in  a closed-loop  phase-compensation 
system  Is  small,  the  spatial  correlation  of  the  aperture  field  may  be 
written 


(r  )u  (r 
an  a an  a 


(4-24) 


Substituting  Eq  (4-24)  Into  the  result  for  the  Interferometer  output  signal 
derived  in  Cnapter  III  (Ref  Eq  (3-30)),  the  x-channel  output  signal  may  be 
calculated  and  Is  written  on  the  next  page  In  Eq  (4-25) . Henceforth  It 
will  be  assumed  that  conditions  are  such  that  the  phase  which  Is  measured 

by  the  Interferometer  represents  the  wavelength- Independent  phase  <Ji_  (r  ) . 

Tq  a 

Eq  (4-25)  gives  the  x-channel  output  signal  of  the  Interferometer 
wavefront  sensor  due  to  a broadband.  Incoherent  source  viewed  through 
atmospheric  turbulence.  Each  temporal  mode  enclosed  by  the  first  raal 
operator  Is  a result  of  the  coherent  addition  of  the  zero  and  +1  dif- 
fracted field  orders,  and  represents  the  positive-shear  term  for  each 
temporal  mode.  The  sum  of  these  terms  expresses  the  superposition  of  all 
positive  shear  modes  for  the  white-light  interference  pattern.  Similarly, 
each  mode  enclosed  by  the  second  real  operator  Is  the  result  of  the  zero 
and  -1  field  orders  and  Is  called  the  negative-shear  mode. 

A few  general  comments  can  be  made  about  the  wavefront  sensor  output 
signal  at  this  time: 
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(1)  All  modes  of  the  detector  signal  are  modulated  at  the 


fundamental  frequency  f^  . If  the  wideband  temporal  process  w(t)  Is 
characterized  over  an  Interval  T much  shorter  than  the  turbulence- 
induced  phase  fluctuations,  the  phase-front  across  the  aperture  may  be 
considered  frozen,  and  no  additional  frequency  modulation  Is  Introduced. 

(2)  For  each  positive-shear  term,  the  total  spatial  effect  of 
the  source  radiance  Is  confined  to  an  Integral  over  the  source  Intensity 
distribution.  For  the  negative-shear  terms,  the  effect  Is  represented 
by  a similar  Integral.  The  Integrals  are  Independent  of  temporal  mode, 
so  that  the  spatial  chracterlstlcs  of  the  source  are  entirely  separate 
from  the  broadband  temporal  effects  of  the  source.  This  result  Is  due 
to  the  coherence  separability  of  the  source  field.  The  detailed  effects 
of  the  source's  spatial  content  are  discussed  at  length  later  In  this 
chapter. 

(3)  For  positive  and  negative-shear  modes,  the  total  effect 
of  the  broadband  nature  of  the  field  Is  contained  In  a sum  over  all  con- 
tributing temporal  modes.  Each  term  In  the  sum  Is  weighted  by  the  power 
spectrum  of  the  broadband  process  w(t)  . 

In  the  next  section,  the  effects  of  the  complex  source  radiance  and 
the  broadband  spectrum  are  separated,  and  the  wavefront  sensor's  phase 
measurement  Is  derived  from  the  Interferometer's  detector  signal  expressed 
In  Eq  (4-25) . 

Derivation  of  the  Wavefront  Sensor  Phase  Measurement 

The  interferometer  wavefront  sensor's  phase  measurement  Is  contained 
in  the  detector  signal  <c^(r^,t)  expressed  In  (4-25).  To  make  the  ana- 
lysis more  Intuitive,  the  spatial  and  broadband  effects  of  the  source  can 


i 

*■ 
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be  separated  and  written  In  more  compact  functional  forms.  Following 
this  separation,  the  derivation  of  the  Interferometer  phase  measurement 
follows  easily. 

Contribution  Due  to  Source  Radiance  Distributions.  As  pointed  out 
In  the  previous  section,  the  total  effect  of  the  source's  spatial  distri- 
bution Is  contained  In  two  similar  Integrals  over  the  source  radiance 
(Intensity)  distribution.  Define  J^(f^,fy)  to  be  the  spatial  Fourier 
transform  of  the  source  radiance  distribution 


J (f  ,f  ) - F {I  (x  ,y  )} 
s X y X y s s ■'s 

^ S'^S 


(4-26) 


Note  that  the  Integral  over  the  source  distribution  In  the  negat Ive-shear 
modes  can  be  written 


- f 


dx 


- F {I  (x  ,y  )} 

Vs  8 8 3 

where  the  Fourier  transform  Is  evaluated  at  spatial  frequencies 


(4-27) 

(4-28) 


f “ ^ ■ 0 . 

X XqZ  y 

Let  the  transform  of  the  source  radiance  distribution  be  written 
In  polar  form  so  that 


where 


J (f  ,f  ) - |j  (f  ,f  ) 
s * s X*  y^ 


(4-29) 

(4-30) 


la  the  magnitude  and 


i>  (f  ,f  ) " arg[j  (f  ,f  )] 
'^s  x’  y'  s x'  y •* 


(4-31) 
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Is  the  phase  of  the  transform  of  the  source  radiance  distribution.  Using 
Eq  (4-31) , the  Integral  (4-27)  may  be  represented 


(4-32) 


where  f 

X 


Msn 
XqZ  * 


Since  the  intensity  distribution  of  the  source  I (r  ) is  a real 

s s 

function  of  the  source  coordinates,  the  inverse  Fourier  transform  of 
the  source  distribution  is  equal  to  the  complex  conjugate  of  the  Fourier 
transform  of  the  source  distribution  (Ref  31:28): 


^X-V  (‘-w 

* y 

Using  this  notation,  the  Integral  over  the  positive-shear  modes  in 
Eq  (4-25)  may  be  written 


3;(f.,o)  - 


Js  ^ ^X  ’ '‘'s  ^ ^x  ’ ^ 


(4-34) 

(4-35) 


Contribution  Due  to  Broadband  Emission.  Consider  the  sum  over  all 


the  positive-shear  temporal  modes  defined  below: 


Q(M(r^-so))  C^(n)S„(7)exp[jir  ] 


w'  T 


« T ^ T 
for  - Y < t 


expCj  (M(r^-SoC ‘^n)  ) ) ] 


(4-36) 
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Using  the  blnonlal  expansion  for  C (n)  and  approximating  the  sum  by  the 
first  two  terms  allows  (4-36)  to  be  written 


Q(M(r^-So))  - exp[j  i I 


exp[j27r(-^)  exp{j4.  [M(r  -sod  + f\))]} 


To‘ 


foT' 


for  (4-37) 

Performing  the  same  approximation  on  the  negative-shear  modes  yields 

Q(M(r^+8o))  » exp[j  -^]  i I c2(n)S^(^) 

n 

( _ _ 

exp[j2ir(-^)^]  exp{j^^^[M(T^+Sod+-^))]  } 


- T ^ T 

for  - y < t <-2 


(4-38) 


Furthermore,  let  the  positive-  and  negative-shear  terms  be  written  in 
polar  form  so  that 

_ _ M^s^ 

Q(M(rd-So))  - e*PCj’f"X^]^(M(r^-so)) 


exp[j'l'^(M(r^-8o))] 


(4-39) 


( 
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and 
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Q(M(r^+so))  - exp[jir-jp£-]  A^(M(r^+So)) 


exp[j((i^(M(F^+8o))] 

(4-40) 

It  Is  shown  In  a later  section  that  for  white-light.  Incoherent  source  fields, 
the  lateral  shear  AC  Interferometer  causes  the  diffracted  aperture  field 
•orders  to  be  spatially  filtered.  The  amplitude  and  phase  of  the  filtered 
aperture  field  are  Identified  by  "hats." 

Interferometer  Phase  Measurement.  Substituting  the  defined  relations 
(4-32) , (4-35) , (4-39) , and  (4-40)  Into  the  x-channel  output  signal  given 
in  (4-25)  yields 

exp[-j2Tr*^  (*j“‘^)]exp[-jTr^(f^,0)] 
exp[j(|i^(M(r^-Io))-j'l'^^(Mr^)]  } 

TTXqZ)^  Jg(f^,0)A^(M(7^+Io))  Re{exp[j(2iTf^t-|)] 

2 

exp[+J2iT^^  (Xj-^)]exp[+ji|)^(f^,0)] 

exp[  j ()>^(M(r^+8c)  ) -1  (to^)  ] } 

(4-41) 
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If  the  amplitudes  of  the  spatially  filtered  orders  are  assumed  to  be 
equal  over  the  distance  of  shear  for  each  point  in  the  detector  plane 

A^(m7^)  = A^(M(7^-Io))  (4-42) 

= ^(M(7^+¥o))  (4-43) 

both  positive- and  negative-shear  terms  In  Eq  (4-41)  may  be  combined: 


>td(rd,t) 


2m2 


■ir(AoZ)  cos[i(.>;+^‘)-<^^^-0j 


sln[2Trf^t  +j  (<|i^  - ({i^) 


where 


- ♦^(M(r^-so)) 

" (|>^(M(r^+so)) 

♦to  ■ ♦i„<“d> 

9,  - ^ 

«d  * ¥ *‘^“»*d 


(4-44) 

(4-45) 

(4-46) 

(4-47) 

(4-48) 

(4-49) 


If  the  amplitudes  of  the  spatially  filtered  orders  are  not  equal  over 
the  distance  of  shear,  the  resulting  detector  signal  Is  represented  as  the 
sum  of  two  phasors  as  discussed  In  Chapter  III  (see  Figure  8) . In  Eqs 
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(4-42)  and  (4-43),  however,  the  filtered  amplitudes  A^(M(r^-8o)) 
and  A^(M(r^+so))  represent  the  amplitude  of  the  diffracted  and  fil- 
tered aperture  phase  exp[j({i_  (r  )]  Introduced  by  atmospheric  turbulence. 

Tq  a 

Since  this  phase  term  Is  assumed  to  have  unit  modulus  for  all  points  In 
the  aperture  plane,  any  sheared  and  filtered  version  of  this  filter  will 
be  equal  to  any  other  filtered  version  so  long  as  the  filtering  effects 
are  Identical.  In  a later  section  It  Is  shown  that  the  amplitude  of  the 
filtered  aperture  field  Is  In  fact  the  same  for  either  the  m - 1 or 
m ■ -1  diffracted  order. 


Therefore,  Eq  (3-44)  represents  the  Interferometer  output  signal  for 
x-shear  for  spatially  Incoherent,  white-light  sources  viewed  through  atmo- 
spheric turbulence  subject  to  the  operating  conditions  discussed  earlier. 
Note  that  all  effects  of  the  source's  radiance  distribution  are  contained 
In  a multiplicative  attenuation  term  J (f  ,0)  and  an  additive  phase 

8 X 

term  i|;  (f  ,0)  . The  effect  of  the  diffraction  of  broadband  light  affects 

8 X 

the  phase  measurement  by  filtering  the  measured  phases  and  (|i~  , 

and  attenuating  the  carrier  by  an  amount  A^(Mr^)  . The  effects  of  source 
radiance  distribution  and  broadband  emission  spectrum  are  discussed  In 
the  next  two  sections. 

As  defined  In  Eq  (4-31),  <|'(£^«0)  Is  the  phase  of  the  spatial  Fourier 

transform  of  the  source  radiance  distribution  evaluated  at  spatial  fre- 
quencies f * and  f >'  0 . Since  the  source  distribution  Is  always 

X AqZ  y 

positive  and  real,  i|*^(f^,0)  Is  an  odd  function  of  f^  and  represents 

a constant  phase  offset  In  Eq  (4-44)  for  any  non-zero  shear  value.  The 

size  of  this  phase  error  depends  on  the  specific  radiance  distribution 

of  the  source  and  the  value  of  f ■ . 

X \qZ  • 
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In  addition  to  the  effects  mentioned  above,  two  additional  phase  teras 
0^  and  6^  appear.  Whereas  , and  are  representative  of 

the  turbulence-induced  phase-front  Introduced  at  the  aperture  plane,  the 
additional  phase  terms  0^  and  0^  represent  the  quadratic  phase  of  a 
spherical  wave  emlnatlng  from  a point  source  located  at  r^  ~ (0,0)  in  the 
source  plane. 

Measurement  of  Turbulence-Induced  Phase.  The  total  phase  of  the 

sinusoidal  signal  modulated  at  frequency  f^  Is  the  sum  of  three  phases: 

(1)  the  differential  aperture  phase  due  to  atmospheric  turbulence 

"2  > (2)  the  phase  of  spatial  Fourier  transform  of  the  source  radiance 

distribution  ij;  (f  ,0)  , and  (3)  a quadratic  phase  0,  . If  the  interfero- 
s X a 

meter  is  set  up  to  measure  the  slope  of  the  turbulence-induced  aperture 
phase,  the  post-detection  processing  of  the  Interferometer  output  signal 
must  be  designed  to  extract  the  phase  due  to  the  first  effect  ((^^-(^^}  , 
alone.  Assuming  signal  processing  has  extracted  the  differential  phase 

- Ij  (A-50) 
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visibility  Limitations  of  Extended  Source  Distributions 

In  Eq  (4-44)  it  Is  seen  that  the  total  effect  of  the  source  radiance 
(Intensity)  distribution  on  the  aperture  phase  measurement  Is  contained 
In  the  amplitude  and  phase  of  J (f  ,0)  defined  by  Eq  (4-27): 

8 X 

*•  r 

Jg(f^.O)  - jd7gIg(7g)expC-j2irfgXj  (4-53) 


where  f * . 

X XqZ 

The  Integral  over  x^  represents  the  spatial  Fourier  transform  of  the 
x-variatlons  of  the  source  Intensity,  evaluated  at  a spatial  frequency 
Integral  over  y gives  the  total  Intensity  along  the 

X ^0^  ® 

y-directlon  for  each  x 

s 

The  result  given  In  Eq  (4-53)  represents  the  mutual  Intensity 
(Ref  14:508)  of  light  at  two  points  In  the  aperture  plane  due  to  a 
quasl-monochromatlc  source  emitting  radiation  near  wavelength  Xg  . The 
points  In  the  aperture  are  separated  In  the  x-dlrectlon  by  distance 
Ax^  ■ Msg  , with  no  separation  in  the  y-dlrection. 

The  magnitude  of  J^(f^,0)  corresponds  to  the  degree  of  coherence 
of  the  aperture  field  (Ref  14:510)  which  Is  defined  to  be  the  absolute 
value  of  the  normalized  Fourier  transform  of  the  source  Intensity  dis- 


tribution: 


vjr.r:  ) - 


/drj„(rjexp  - j <Vs'^ays> 


a a a 


fdr  I (r  ) 

^ BBS 


(4-54) 


where  P (r  ,t')  Is  the  degree  of  coherence  of  the  light  at  two  points 
AAA 

r^  and  r'  In  the  aperture  plane. 


V 


1 
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If  the  amplitudes  of  the  two  sheared  field  orders  are  equal — as 


expressed  In  Eqs  (4-42)  and  (4-43) — then  J^(f^,0)  Is  related  to  the 
fringe  visibility  (Ref  14:507-508)  of  the  Interference  fringes  at  the 
detector  plane: 


y(Mso) 


J (f  ,0) 

s X 

Jdr  I (r  ) 
■'  s s s 


(4-55) 


where  y(Mso)  Is  the  fringe  visibility  of  the  Interference  fringes  for 
two  points  In  the  aperture  plane  separated  by  distance  Msg  In  the 
x-dlrectlon  with  no  separation  In  the  y-dlrectlon,  and  0 _<  yCMsg)  _<  1 
for  any  shear.  Note  that  for  the  case  of  a point  source,  the  fringe 
visibility  Is  unity  for  all  shear  values.  For  any  extended  source, 
however,  the  fringe  visibility  depends  on  the  specific  source  radiance 
distribution  and  the  shear  distance  chosen. 

Clearly,  J^(f^,0)  Is  proportional  to  the  amplitude  of  the 
Interferometer  output  signal  for  x-shear,  and  must  be  maintained  at  a 
significant  value  for  the  Interferometer  signal  to  be  detected  and  the 
phase  Information  to  be  extracted.  The  effects  of  shear  distance  sg 
on  the  fringe  visibility  are  considered  below  for  the  case  of  uniform  and 
complex  extended  source  radiance  distributions. 

Uniform  Source  Distributions.  The  visibility  of  the  Interferometer 
fringe  pattern  for  x-shear  Is  the  normalized  Fourier  transform  of  the 
source  Intensity  distribution  as  expressed  In  Eq  (4-55).  As  an  example, 
consider  uniform  source  with  a square  Intensity  distribution  with  width 
W meters  and  Intensity  Ig  watts  per  meter  squared.  The  Interferometer 
fringe  visibility  can  be  calculated  from  Eq  (4-55) : 
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4. 


(4-56) 


Y^(Mso)  - 

where  sinc(x)  * sin(x)/x  . The  fringe  visibility  of  the  square  source 

is  plotted  as  a function  of  shear  distance  Msq  in  Figure  12. 

For  a uniform  square  source,  then,  the  interference  fringe  visibility 
generally  drops  off  as  the  shear  distance  Increases.  To  guarantee  that 
significant  fringe  visibility  is  maintained,  the  shear  is  selected  so  that 
Yj^(Msq)  is  evaluated  within  the  main  lobe  of  the  sine  (•)  function: 

Mso  < ^ (4-57) 

where  a >•  W/Z  is  the  angular  subtense  of  the  source. 

For  uniform,  circular  source  distributions  the  fringe  visibility  is 
related  to  a first  order  Bessel  function,  and  a result  similar  to  Eq  (4-57) 
for  circular  sources  can  be  shown  to  be 

Msn  < 1.22  (4-58) 

^ a 

where  a ~ D/Z  is  the  angular  subtense  of  the  circular  source. 

Therefore,  the  restriction  placed  on  shear  distance  for  circular 
sources  with  angular  subtense  a * D/Z  is  on  the  same  order  as  (4-57) 
for  a square  target  of  equal  angular  subtense. 

Complex  Source  Distributions.  Sources  of  interest  to  imaging  systems 
generally  contain  some  degree  of  fine  structure.  Sources  with  structure 
contain  more  high  spatial  frequency  information  than  uniform  sources, 
and  cause  the  fringe  visibility  function  y(Mso)  to  have  significant  value 

( 

i 

85 

f 

( 


sinc[iT(  Y^)w] 
AnZ 


4 


AO-A06H  HOH  AIR  FORCE  INST  OF  TECH  KRIGHT-PATTERSON  AFB  OHIO  SCH— ETC  F/S  20/6 

COHERENCE  PROPERTIES  OF  BROADBAND  OPTICAL  FIELDS  WITH  APPLICATI— ETC (IM 
OeC  TB  PS  IDCU. 

UNCLASSIFIED  AFIT/BC0/Ce/7t-l  NL 


f 


for  large  shear  values.  Although  the  visibility  effects  of  each  complex 
source  must  be  considered  separately,  the  general  effect  can  be  Illustrated 
by  the  following  example. 

Consider  a square  source  of  width  W having  a periodic  Intensity 
distribution  similar  to  a bar  pattern  with  variations  along  the  x-dlrectlon 
as  shown  In  Figure  13.  Note  that  the  visibility  function  YjCMso)  contains 
weighted  copies  of  the  visibility  function  of  a square  source  VyCMsg) 
located  at  sq  ■ m^  , where  m Is  an  integer.  It  Is  seen  that  periodic 
source  distributions  have  a "coherence  modulation"  affection,  due  to  the 
spatial  frequency  modulation  properties  of  periodic  source  functions;  and 
the  fringe  visibility  can  have  significant  value  for  shear  distances  much 


larger  than  those  allowed  by  uniform  sources. 

<■ 

J , 

F 
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Figure  13.  (a)  Periodic  Source  Dist 

Fringe  Visibility  Function  yCMsq)  f' 


since  the  total  radiant  power  for  any  source  Is  finite,  the  power  Is 
distributed  among  the  many  lobes  of  the  fringe  visibility  function,  and 
the  magnitude  of  the  fringe  visibility  at  any  shear  may  be  negllblble. 

In  any  practical  sense,  then,  the  measurable  fringe  visibility  Is  greatly 
dependent  on  total  source  radiance  and  the  complexity  of  the  source  dis- 
tribution . 


Special  Filtering  Effects  of  Broadband  Light 

Consider  the  approximation  to  the  sum  of  the  posit Ive-shear  temporal 
modes  written  previously  In  Eq  (4-37): 


( 


Q(M(rj-so)) 


exp[je.]  7 I C^(n)S^(Y)exp[j2iT( 


2.2 


M-^s 

Ic 


■)f] 


p{j*~  [M(r.-ao(l+7^))]} 


'To 


foT 


- T ^ T 
for  " "2 ^ ^ Y * 


(4-59) 


Taking  the  spatial  Fourier  transform  of  (4-59)  with  respect  to  the  aper- 
ture plane  coordinates  x - Mx.  and  y - My.  yields 

A Q A u 

y {Q(M(F^-8o))}  - exp[je^]*^^(f^,fy)exp[-j2irf^M8o] 


r 


T T 
for 


1 ^ ^0 
Y I C^(n)S^(  Y)cxp[-J2Tr(  — 
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(4-60) 


where 


>_  (f  ,f  ) ■ F {exp[j<>_  (r  )]} 
To  x’  Y X y/  a"-' 

& O 


(4-61) 


For  broadband  power  spectra  and  long  process  characterization  Intervals 
the  sum  over  temporal  modes  may  be  approximated  by  an  Integral: 


^x  y tQWr^-So))}  - exp[j9g]4>^^(f^,fy)exp[-j2Trf^Mso] 


jdv  C^(v)S^(v)exp[-j2Trc^v] 


(4-62) 


where 


r . Ikfl.  ff  _ \ 

'^+  fo  '■  X XnZ^ 


(4-63) 


and  V has  been  substituted  for  — as  the  Integration  variable.  The 
Integral  over  v In  Eq  (4-62)  represents  the  first-order  spatial  filter 
for  the  m ■ +1  diffracted  aperture  field  order  given  by  Eq  (3-87) 
where  the  complex  degree  of  coherence  of  white-light  aperture  field  has 
been  calculated  explicitly. 

Following  Eq  (3-88)  the  first-order  spatial  filter  affecting  the 

m " +1  diffracted  aperture  phase  exp[j^  (r  )]  Is 

Jo  a 


HJC.)  - F^{52(v)s„(v)} 

▼ TV  W 


(4-64) 


For  typical  compensated  Imaging  applications,  sources  are  distant 
compared  to  the  shear  distance  Mso  . Assuming  Is  small  compared 

to  the  spatial  bandwidth  of  H 


®H  x5z 


(4-65) 


the 


where  Is  the  effective  bandwidth  of  the  spatial  filter  , 

first-order  spatial  filter  can  be  written 

f^Mso 

irhere  Z “ — as  defined  In  Eq  (3-60) . Using  results  from  Chapter  II 

(Ref  Eq  (2-71)),  Eq  (4-66)  can  be  written  In  terms  of  the  first  two  deri- 
vatives of  the  temporal  correlation  function  R^(t)  : 

V')  ^ ^ KM  - ^ KM  (‘-6') 


Similarly,  the  sum  over  the  negative-shear  temporal  modes  can  be 
approximated  by  an  Integral,  and  a filtering  function  for  the  m > -1 
diffracted  field  order  can  be  defined.  Following  Eqs  (3-90)  and  (2-70), 


H (C)  - F“^{e2(v)S  (v)} 
• V w 


(4-68) 


irfo  3? 


1_ 

(2irfo)'2‘ 


- v«> 


(4-69) 

(4-70) 


Therefore,  Eq  (4-67)  and  Its  conjugate  (4-89)  represent,  at  least  to 
first  order,  the  form  of  the  spatial  filtering  Introduced  Into  the  Inter- 
ferometer phase  measurement  due  to  the  diffraction  of  broadband  light  for 
the  m ■ -t-l  and  m ■ -1  diffracted  orders.  The  spatial  filtering  func- 
tion Is  similar  to  the  "windowing"  function  disctissed  In  Chapter  II 

with  regard  to  coherence  length  (Ref  Eq  (2-70) ) . 

Indeed,  If  the  power  spectrum  S^(f)  sufficiently  narrow  so  that 
C^(n)  may  be  closely  approximated  by  unity,  the  filtering  functions 
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and  H take  on  the  form  of  the  temporal  correlation  function: 


H^(0  - H_(C)  - R^(0  (4-71) 

Conalder,  for  example,  a bandllmlted  power  spectrum  S.  (f)  discussed  In 

D 

Chapter  III  (Ref  Eq  (3-70)).  The  effective  bandwidth  of  the  resulting 
spatial  filter  Is  on  the  order  of 


(111) 

MSg  ^Af^ 


(4-72) 


as  defined  In  Eq  (3-71) . If  the  bandwidth  of  the  spatial  Fourier  transform 
of  the  turbulence-induced  aperture  phase  4^(f^,f^)  Is  much  less  than 


B,  « B 


H 


(li) 

Msq  ^Af^ 


(4-73) 


where  B Is  the  spatial  bandwidth  of  either  real  or  Imaginary  quadrature 

of  exp£j^_(r  )]  , then  the  spatial  filters  Introduce  negligible  filtering 
Ig  a 

of  the  aperture  phase  field  orders  contributing  to  the  AC  Interferometer 
pattern,  and  cause  little  distortion  In  the  measured  aperture  phase. 

In  summary,  then,  for  truly  broadband  optical  fields  with  arbitrary 
power  spectra,  the  form  of  the  spatial  filtering  must  be  calculated  from 
tha  discrete  sum  over  temporal  modes  as  In  Eq  (4-36).  To  first  order, 
however,  the  spatial  filtering  can  be  calculated  from  derivatives  of  the 
broadband  field  correlation  function  as  In  Eqs  (4-67)  and  (4-69) . 
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V.  CCWCLUSION 


( 


Summary  of  Results 

Coherence  Model  for  Broadband  Fields.  A free-space  propagation 
model  for  broadband  optical  fields  was  developed  based  on  a Karhunen- 
Lodve  (KL)  expansion  of  the  time-varying  portion  of  a coherence  separable, 
broadband  optical  envelope.  Due  to  the  linearity  of  the  free-space  chan- 
nel, each  temporal  mode  of  the  field  expansion  could  be  propagated  Indi- 
vidually, and  the  output  broadband  field  could  be  computed  as  simply  the 
superposition  of  all  the  Individually  propagated  field  modes.  It  was 
found  that  for  broadband  optical  fields  characterized  over  long  time 
Intervals,  the  eigenfunctions  of  the  KL  expansion  were  approximated  by 
the  complex  exponentials  of  a Fourier  series  expansion  over  the  same  time 
Interval.  The  eigenvalues  of  the  modal  expansion  were  found  to  be  samples 
of  the  temporal  power  spectrum,  sampled  at  the  harmonic  frequencies  of 
the  Fourier  series  expansion.  Since  the  expansion  coefficients  of  a KL 
expansion  are  uncorrelated,  the  calculation  of  output  field  correlation 
was  facilitated.  The  output  field  correlation  for  coherence  separable 
source  fields  was  stated  for  special  cases  of  source  field  coherence 
In  Table  I.  The  output  field  correlation  due  to  coherence  separable 
source  fields  was  found  to  be.  In  general,  not  coherence  separable 
(Ref  Eq  (2-47)). 

Analysis  of  White-Light  Interferometer.  Application  of  the  broad- 
band propagation  model  to  the  lateral  shear  heterodyne  Interferometer 
showed  that  the  position  of  the  Interference  fringes  detected  by  the 
Interferometer  was  dependent  on  wavelength.  This  effect,  caused  by  the 
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diffraction  of  the  white-light  aperture  field  by  the  interferometer’s 
rotating  grating,  was  shown  to  spatially  filter  the  m « 1 and  m - -1 
diffracted  field  orders  contributing  the  Interference  pattern.  To 
first-order,  the  form  of  the  spatial  filter  was  shown  to  be  given  by  a 
Fourier  transform  of  the  product  of  the  broadband  temporal  power  spectrum 
and  a shifted  version  of  the  complex  degree  of  coherence  of  the  aperture 
field  (Ref  Eqs  (3-87)  and  (3-89)).  The  measurement  made  by  the  inter- 
ferometer was  found  to  be  related  to  the  phase  of  the  spatially  filtered 
aperture  field  envelope.  The  Interferometer  wavefront  measurement  was 
specialized  for  spatially  coherent  aperture  fields  for  large  and  small 
shear  (Ref  Eqs  (3-83)  and  (3-77)). 

Measurement  of  Turbulence-Induced  Phase.  The  broadband  field  coherence 
model  was  applied  to  the  application  of  the  shearing  Interferometer  as  a 
wavefront  sensor  in  a phase-compensated  imaging  system.  It  was  argued 
that  if  the  turbulence-induced  phase  fluctuations  '^xn^^a^  were  correlated 
over  the  wavelength  range  of  the  source  emission  spectrum  or  the  sensor 
operated  in  a closed-loop  mode  where  the  residual  phase  tracking  errors  were 
small,  the  turbulence- Induced  phase  fluctuations  could  be  modeled  as  a unit- 
modulus  phase  screen  for  all  wavelengths  of  interest  (Ref  Eq  (4-24) ) . 

If  the  source  was  immobile  and  sufficiently  diffuse,  the  source  field 
could  be  modeled  as  being  coherence  separable.  This  model  allows  the 
visibility  effects  due  to  extended  source  radiance  distribution  and  the 
spatial  filtering  effects  due  to  broadband  emission  to  be  separated  in 
the  analysis  of  the  interferometer.  The  visibility  of  the  Interference 
fringes  was  found  to  be  proportional  to  the  magnitude  of  the  Fourier  trans- 
form of  the  source  radiance  distribution  as  prescribed  by  the  Van  Clttert- 
Zemieke  theorem  (Ref  14:510)  for  quasi-monochromatlc  sources  (Ref  Eq  (4-28)). 


It  was  found  that  the  phase  measured  by  the  Interferometer  for  this 
application  was  due  to  the  sum  of  three  effects  (Ref  Eq  (4-A4)):  (1) 

the  turbulence-induced  aperture  phase  perturbation,  (2)  the  phase  of 
the  spatial  Fourier  transform  of  the  source  radiance  distribution,  and 
(3)  the  quadratic  phase  due  to  a point  source  at  the  source  plane  mea- 
sured In  the  aperture  plane.  In  order  to  measure  the  aperture  phase 
due  to  atmospheric  turbulence,  signal  processing  must  be  designed  to 
extract  the  first  phase  effect  from  the  Interferometer  output  signal. 

For  typical  compensated  Imaging  situations  It  was  found  that  to  first 
order,  the  spatial  filtering  could  be  expressed  In  terms  of  the  first  two 
derivatives  of  the  broadband  temporal  correlation  function  (Ref  Eqs  (4-67) 
and  (4-69)). 

Discussion  and  Suggestions  for  Further  Study 

In  Chapter  II  a coherence  model  for  the  propagation  of  broadband 
optical  fields  was  developed  from  a modal  expansion  of  the  time-varying 
part  of  the  complex  field  envelope.  Although  the  model  was  specifically 
applied  to  the  analysis  fo  a white-light  shearing  Interferometer  In  this 
paper,  the  application  of  this  model  may  be  extended  to  many  broadband 
optical  or  Infra-red  systems  for  which  the  Input  fields  can  be  assumed 
to  be  coherence  separable.  For  exaiqtle,  the  free-space  propagation 
model  (Eq  (4-42)),  the  apochromatlc  lens  result  (Eq  (3-7)),  and  the 
rotating  diffraction  grating  model  (Eq  (3-13))  are  particularly  appli- 
cable to  optical  Interferometric  and  measurement  systems.  Indeed,  a 
broadband  optical  field  can  still  be  represented  as  a modal  expansion 
If  It  la  due  to  a source  which  can  be  assumed  to  be  coherence  separable 
(Ref  Eq  (2-42))  and  the  effects  of  the  propagation  medium  can  be  considered 
to  be  frozen  In  time  (Ref  Eq  (4-8)). 
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In  Chapter  III  the  lateral  shear  AC  Interferometer  was  analyzed 
for  broadband  aperture  fields  with  arbitrary  spatial  coherence.  Several 
assumptions  were  made  which  may  further  qualify  the  application  or  degrade 
the  predicted  performance  of  the  wavefront  sensor.  First  of  all,  all 
diffraction  effects  due  to  finite  lens  apertures  and  field  stops  were 
Ignored.  All  optics  were  assumed  to  be  perfectly  transmissive  for  all 
wavelengths  of  Interest  and  the  lenses  were  assumed  to  be  aberration-free 
and  perfectly  apochromatlc.  These  effects  must  be  studied  In  detail, 
or  at  least  considered,  for  any  practical  application  of  a physically 
realizable  wavefront  sensor. 

The  Interferometer  output  signal  ^^(r^,t)  calculated  In  Chapters 
III  and  IV  (Ref  Eqs  (3-43)  and  (4-44))  represents  the  time-varying 
Interference  pattern  at  each  point  In  the  detector  plane  modulated  at 
f^  Hertz.  No  model  for  the  photon/electron  conversion  of  this  white- 
light  Intensity  pattern  Into  an  electrical  signal  was  proposed,  nor  was 
a signal  processing  scheme  suggested  for  the  extraction  of  the  aperture 
phase  measurement  A(P  from  the  detector  plane  output  signals.  It  Is 
recommended  that  a model  for  the  detection  of  broadband  radiation  be 
Incorporated  Into  the  analysis  of  the  Interference  pattern  (Eq  3-34) 
to  determine  what  additional.  If  any,  filtering  effects  are  Introduced 
Into  the  wavefront  measurement.  Once  these  detected  Interferometer 
output  signals  are  determined,  several  electronic  processing  schemes 
can  be  considered  to  provide  the  best  aperture  phase  measurement 
(e.g..  Refs  32,  33,  and  34). 

It  was  found  that  the  shearing  Interferometer's  operation  for 
broadband  fields  causes  the  m ■ +1  and  m -1  diffracted  aperture 
field  orders  contributing  the  detector  plane  Interference  pattern  to 
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be  spatially  filtered.  Whereas  the  filtering  affects  both  real  and 
Imaginary  quadratures  of  the  complex  aperture  field  envelope,  the  phase 
measurement,  derived  from  the  Interferometer  output  signal  (Eq  (3-47)), 
represents  the  phase  of  a filtered  complex  field.  Except  for  the 
special  filtering  cases  studied — extremely  wide  filter  (Eq  (3-73)  and 
extremely  narrow  filter  (Eq  (3-83)) — the  effect  of  spatial  filtering 
on  the  detected  aperture  phase  Is  extremely  difficult  to  determine.  The 
effects  of  filtering  on  the  measured  phase  of  phase  modulated  systems  have 
been  studied  (e.g..  Refs  35,  36,  and  37),  but  are  beyond  the  scope  of  this 
paper. 

For  the  analysis  of  shearing  Interferometer  wavefront  measurement 
It  was  assumed  that  either  the  phase  of  the  complex  aperture  field 
envelope  was  the  same  for  all  wavelengths  of  light  or  the  Interferometer 
was  working  In  closed-loop  with  a real-time  phase-correction  system. 

The  end  result  was  that  the  Interferometer  adequately  measured  a phase 
which  Is  common  for  all  wavelengths  when  the  residual  phase  errors  are 
small.  These  assu]iq>tlons  were  motivated  by  the  fact  that  for  the  compen- 
sated Imaging  system  considered  in  Chapter  IV  (Ref  7) , atmospheric  phase 
correction  Is  accomplished  by  a single  corrector — thus  Introducing  the 
sane  compensation  for  all  wavelengths.  For  more  general  active  optic 
applications  In  which  the  wavefront  correction  is  wavelength-dependent 
(e.g..  Ref  12),  the  decomposition  of  the  aperture  phase  into  wavelength- 
dependent  temporal  modes,  and  the  effect  of  this  decomposition  on  the 
Interferometer  measurement,  may  be  studied. 
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